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Preface 


This  work  was  performed  by  Mr.  David  L.  Thirkill  while  attending  the  University  of 
California  at  Davis.  The  Hydrologic  Engineering  Center  (HEC)  contracted  with  the 
University  for  this  research  and  this  report  was  Mr.  Thirkill’s  M.S.  Thesis.  Mr.  Arlen 
Feldman  of  the  HEC  was  on  Mr.  Thirkill’s  thesis  committee. 

The  ADAPT  software  used  in  this  research  is  no  longer  available.  It  was  proprietary 
software  of  W.E.  Gates  and  Associates.  The  Corps  has  not  continued  with  its  usage  because 
of  the  many  new  commercial  ^tems  now  available. 

The  modelling  concepts  discussed  herein  with  the  ADAPT  software  could  be 
duplicated  with  currently  available  vector-based  geographic  information  interns  (GIS)  e.g. 
ARC/INFO  TIN  model.  HEC  is  continuing  research  in  the  area  of  GIS-based  hydrologic 
analysis. 
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1.  IHTRODUCTIOM 


1.1  Watershed  Modellnc  Using  Geographic  Information  Systems 

The  planning  of  water  resources  projects  relies  heavily  on 
geographic  Information  describing  river  basins.  Information 
about  topography,  land  \ise.  vegetative  cover,  soil  type  and 
erodlblllty  are  needed  In  rainfall -nmoff  modeling,  flood  damage 
determination,  soil  erosion  studies  and  water  quality  studies.  It 
is  often  necessary  to  work  manually  with  this  data  to  derive 
input  for  various  simulation  models  and  resource  studies. 
Geographic  Information  Systems  (GIS's)  have  been  developed  to  take 
advantage  of  the  data  handling  capability  of  digital  computers 
enabling  more  detailed  modeling  and  easing  the  burden  of  much  of 
this  hand  work. 

A  basic  problem  confronting  water  engineers  has  been  how  to 
handle  the  heterogeneity  In  the  geographic  characteristics  within 
a  basin.  Because  of  their  efficiency  in  handling  data.  GIS’s 
have  been  applied  to  this  problem.  A  GIS  allows  engineers  to 
model  the  hydrologic  diversity  within  a  watershed  with  a 
resolution  dependent  only  on  the  size  of  the  elements  chosen.  It 
Is  generally  felt  that  modeling  smaller,  more  homogeneous  areas 
yields  a  more  accurate  simulation  (1,6,20,22,23,25).  With  a  GIS, 
less  averaging  is  required  and  greater  use  of  readily  available 
physical  data  Is  accomplished.  Derivation  of  routing  and  runoff 
coefficients  Is  therefore  based  on  a  more  accurate  physical 
model. 

A  GIS  can  also  provide  the  basis  for  modeling  the  hydrology 


of  ungaged  river  basins  and  for  studying  the  hydrologic  isqiact 
of  physical  changes  (such  as  urbanization)  within  a  river  basin. 

The  conputer  progran  systen,  HECl-ADAPT,  combines  two 
existing  omdels.  ADAPT  Is  a  GIS  that  was  originally  developed  by 
V.E.  Gates  and  Associates  to  aid  In  sewer  design  (34).  HEC*1  is  a 
rainfall -runoff  model  that  was  developed  at  the  Hydrologic 
Engineering  Center  of  the  Corps  of  Engineers  (18) .  These  two 
modeling  systems  are  linked  through  an  Interface  program  called 
HECAD,  also  developed  (under  contract  to  the  HEC)  by  V.E.  Gates 
and  Associates  (34) .  This  report  describes  the  testing  of  the 
HECl -ADAPT  system  for  rainfall -runoff  modeling. 

1.2  Oblectlve.;  of  This  Study 

A  major  objective  of  this  study  Is  to  test  the  ability  of 
the  HECl-ADAPT  system  to  model  rainfall-runoff  processes  on 
ungaged  basins.  Thus,  the  Initial  model  of  each  of  the  two 
basins  studied  is  developed  without  xislng  streamgage  data.  A 
second  objective  Is  to  determine  the  effect  of  model  resolution 
on  the  outflow  hydrograph.  Castro  Valley,  the  smaller  test 
basin.  Is  therefore  modeled  using  two  different  resolutions.  The 
larger  test  basin.  Potter  Valley,  Is  modeled  using  a  single 
resolution.  A  third  objective  Is  to  test  the  flexibility  of  the 
model.  To  evaluate  this,  a  small  urban  basin  and  a  large  non- 
urban  basin  are  iised  In  the  testing  program. 

Before  starting  the  testing  program,  it  was  necessary  to 
accomplish  three  tasks.  The  first  task  was  to  make  sure  all 
necessary  programs  and  hardware  were  available  to  develop  a 


eooplettt  CIS.  The  second  Cask  was  Co  check  each  prograa  Co  make 
sure  ic  worked  properly.  And  Che  chird  cask  was  Co  check 
compuCaClonal  rouCines  Co  be  sure  reasonable  numbers  were  being 
generaced.  This  work  is  described  in  Appendix  A  of  Chis  reporc. 


2.  BACKGROUND 


2.1  Types  of  GIS*« 

Geographic  Information  Systems  (GIS's)  are  data  base  systems 
that  are  used  primarily  for  managing  spatial  geographic  data.  In 
general,  GIS's  have  the  follo%rlng  characteristics:  some  method 
for  entering  and  editing  data  for  the  data  base,  various  systems 
for  displaying  Information  stored  in  the  data  base  and  a 
capability  to  perform  calculations  and  sorting  on  data  in  the 
data  base  (23).  The  types  of  geographic  information  stored  in 
GIS's  are  dependent  on  the  purposes  for  which  the  GIS's  are 
developed. 

There  are  two  basic  t^pes  of  GIS's:  the  polygon  system  and 
the  grid-ceil  system.  The  polygon  system  employs  an  irregular 
polygon  areal  unit  for  spatial  representation.  This  system 
attesipts  to  represent  exact  bovsidaries  of  areas,  points  and 
lines.  Polygon  GIS's  are  used  to  store  maps  in  computers  and  to 
prepare  other  maps  at  different  scales  or  projections  (23). 
Although  the  systems  can  have  high  geographic  fidelity,  they  have 
very  limited  analytic  capablll^  . 

The  grid-cell  system  is  an  Alternative  system  with  much 

» 

improved  analytic  capability  (9,23).  With  the  grid-cell  system, 
the  area  of  interest  is  broken  up  into  square  or  rectangular 
elements  with  various  data  types  and  values  associated  with  each 
element  or  cell.  Analysis  is  usually  done  cell  by  cell. 
Searching,  calculation  of  distance,  production  of  overlay  maps, 
and  stiitabillty  analysis  are  typical  studies  carried  out  vising 


this  Mthodology  (9,23).  Ease  of  manlptilaclon  and  storage  of 
data  using  the  grid  cell  representation  has  resulted  in  wide¬ 
spread  use  of  the  format  as  a  foundation  in  many  GIS's  (9).  A 
major  problem  associated  with  grid-cell  systeu,  however,  is  that 
of  resolution  (7,23,15).  In  order  to  capttire  detail,  a  large 
number  of  grid-cells  are  required,  increasing  both  the  computer 
storage  requirements  and  computing  costs.  Even  with  a  small  grid 
size,  the  system  is  unable  to  precisely  represent  point 
locations,  lines,  or  spatial  boundaries  tislng  the  nodes  of  the 
rectangular  grids. 

A  coid>lnation  system,  the  polygon-to-grid  system,  attenqits 
to  take  advantage  of  the  good  qualities  of  both  the  above 
systems.  Data  is  first  represented  tising  polygons  and  then 
translated  by  computer  to  a  grid  system  for  analysis  purposes 
(23). 

2.2  Digital  Terrain  Models 

Digital  Terrain  Models  (DTM's)  are  considered  to  be  a 
special  type  of  GIS  (23).  In  addition  to  the  usual  attribute 
data  contained  within  the  GIS,  DTM's  also  contain  information  on 
terrain  elevation. 

DTM's  are  normally  produced  using  either  the  rectangular  grid 
or  the  triangular  irregular  network,  a  special  type  of  polygon 
representation  (7).  The  vertices  of  the  triangles  (nodes) 
contain  coordinate -elevation  data  and  the  areas  within  the  lines 
connecting  the  nodes  contain  spatial  data.  The  disadvantage  of 
rectangular  grids  for  terrain  modeling  is  that  the  projection  of 
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the  grid  pattern  onto  a  eooplex  surface  is  warped.  The 
difficulty  of  tislng  the  system  lies  In  tiie  attempt  to  derive 
slopes  and  areas  from  the  warped  quadrilaterals.  Use  of  the 
triangular  Irregular  network  (TIH)  avoids  this  problem  since  the 
vertical  projection  still  yields  a  triangular  plane  that  best 
fits  the  complex  surface  (7) .  Surface  slopes  and  areas  can  be 
computed  easily.  Another  advantage  of  the  TIN  Is  the  ability  to 
vary  the  triangle  size.  Thus,  areas  of  complex  topography  can  be 
accurately  modeled  by  tislng  more  triangles.  Becaxise  of  this 
variable  resolution,  the  TIN  Is  Inherently  more  efficient  than 
the  grid  system  for  modeling  terrain.  The  DTK  being  investigated 
in  this  report  Is  contained  within  ADAPT  (34)  and  is  described 
below. 

2.3  ADAPT 

ADAPT  uses  a  TIN  system  to  store  data.  Terrain  is 
represented  as  a  faceted  surface  with  each  facet  a  triangular 
plane.  Increased  accuracy  of  representation  (resolution)  is 
obtained  where  necessary  by  increasing  the  number  of  triangles. 
This  "variable"  resolution  increases  the  computational 
efficiency  of  the  method. 

Triangle  sides  and  vertices  are  chosen  to  represent 

» 

Important  terrain  features  such  as  ridges,  peaks,  slope  breaks, 
passes  and  streams,  as  well  as  the  natural  boundaries  between 
different  soil  and  land-use  types  and  the  artificial  boundaries 
delineating  political  districts.  Like  the  grid  system,  each 
triangular  element  In  the  network  Is  treated  as  a  homogenoxis 
cell.  Each  cell  contains  Information  such  as  land  use,  soil 


7 


political  Jurisdiction  as  well  as  slope,  slope  direction, 
area  and  elevation.  ADAPT  Inclxides  routines  which  use  Che 
topology  o£  the  DIM  to  deteralne  screen  and  overland  flow 
networks.  The  ADAPT  system  Incorporates  the  good  boundary 
representation  of  the  polygon  systen  while  retaining  the 
analytical  capability  of  Che  grid- cell. 

2.4  Digital  Terrain  Model  Construction  Using  ADAPT 

The  first  step  In  producing  the  DIM  Is  to  delineate  the 
study  area  boundary  on  topographic  maps.  Normally,  7.5  minute 
uses  quads  are  used.  Next,  the  process  of  triangulating  the 
basin  Is  begun  by  overlaying  a  sheet  of  mylar  on  Che  quad. 
Trlaxigles  representing  Che  major  topographic  features  are  drawn 
on  this  overlay.  Each  triangle  shotild  represent  a  xmlform  or 
nearly  uniform  planar  section  of  the  topographic  map.  The 
trlangulatlon  Is  digitized  by  recording  the  coordinates  and 
elevation  of  each  triangle  vertex  using  a  digitizer  connected  to 
a  computer.  A  file  containing  coordinate  and  elevation  data  for 
each  triangle  vertex  ^is  created  this  way. 

This  file  Is  used  Co  build  two  additional  files:  a  triangle 
file  and  a  vertex  file.  Tha  ADAPT  program  used  to  accomplish 
this  Identifies  triangles  wl^  comrK>n  vertices  and  assigns  unique 
vertex  and  triangle  numbers.  It  then  uses  this  Information  to 
determine  triangle  adjacencies  and  to  Identify  basin  boundaries. 
The  resulting  triangle  file  consists  of  vertex  and  adjacency 
information  along  with  coordinate,  elevation,  slc;'>e  and  slope 
angle  data.  Each  vertex  In  the  vezu^x  .lie  contains  a  list  of 


triangl«s  which  share  it  as  a  cobsmxi  vertex  in  addition  to 
coordinate  and  elevation  data.  The  end  restilt  of  this  process  is 
a  DIM.  Soil  and  land  use  data  need  to  be  incorporated  to 
conplete  the  CIS. 

The  additional  data  required  to  produce  the  GIS  are 
typically  derived  from  soil  and  land  \ise  naps  conaonly  available 
from  the  U.S.  Soil  Conservation  Service  (SCS)  and  other  agencies 
involved  in  resource  planning.  This  Information  may  be 
Incorporated  into  the  triangle  file  by  either  digitizing  or  by 
manual  techniques.  Triangles  that  share  the  same  attribute 
characteristics  with  neighbors  (e.g.  same  land  use  or  soils)  are 
aggregated  into  polygons.  The  ADAPT  term  for  these  polygons  is 
"unique  attribute  polygon”  (UAP). 

If  a  digitizing  procedure  is  used,  the  boundaries  of  each 
polygon  are  digitized  and  an  ADAPT  program  assigns  land  use  and 
soils  values  to  individual  triangles  by  determining  what 
triangles  are  internal  to  the  polygon  boundary.  If  done 
manually,  the  Information  is  entered  on  a  triangle-by- triangle 
basis. 

Further  processing  by  the  ADAPT  system  produces  a  network 
file  containing  stream  and  overland  network  data.  This  file 
contains  those  triangle  sides  which  the  program  has  defined  as 
stream  segments  (links).  The  ADAPT  system  automatically  assigns 
stream- link  status  to  a  triangle  side  if  it  is  a  common  side 
between  two  triangles  that  drain  toward  each  other.  Any  other 
triangle  side  can  be  manually  assigned  this  status. 

Coordinate  data,  upstream  and  downstream  vertex  numbers  and 


elevation  date  for  each  link  vertex  are  part  of  the  file.  The 
file  also  contains  channel  roughness  values  and  optional 
Inforaation  describing  stream  cross-sections  (for  normal  depth 
routing),  and  an  overland  network  of  contributing  triangles  for 
each  stream  link. 

2.5  HECAD 

HECAD  Is  the  Interface  program  designed  to  generate  Input 
data  (on  disk)  for  HEC-1  (34)  %islng  Information  stored  within  the 
ADAPT  data  base  and  atixlllary  files.  The  auxiliary  files 
Include:  a  soil  matrix  file  containing  Information  on  each  soil 
type;  a  drainage  network  file  containing  stream  and  overland 
networks  and  channel  Information  required  for  routing;  an 
auxiliary  file  containing  Information  describing  each  channel 
link  for  the  normal-depth  routing  option;  a  sub-watershed 
Identifier  contained  In  the  drainage  network  file;  a  rain  gage 
file  containing  the  ralngage  nuiid>er,  type  and  location;  a 
reservoir  file  containing  routing  characteristics;  a  diversion 
file  describing  location  and  amount  of  diversion;  and  a 
calibration  file. 

The  calibration  file  contains  numerical  values  of 
Infiltration,  roughness  and  percent  imperviousness  as  a  fiinction 
of  land  use  and  soil  hydrologic  group.  HECAD  derives  areally- 
weighted  averages  of  these  parameters  for  each  sub -basin  using 
the  calibration  data  In  the  calibration  file  and  the  soil  and 
land  use  data  in  the  data  base  and  auxiliary  files.  The 
following  sections  describing  HECAD  are  paraphrased  from  the  W.E. 
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6at«a  docuBantaClon  for  HECl -ADAPT  (34). 

2.5.1  Sub-b—ln  Definition 

HEC-1  Is  set  up  to  run  using  the  sub-basin  as  Its  eleo^ntary 
areal  unit.  At  present,  HEC-1  does  not  support  overland  routing 
betveen  sub-baalns.  Therefore,  triangle -to -triangle  routing 
cannot  be  acconpllshed  using  HEC-1.  Soae  amount  of  aggregation 
must  therefore  take  place.  In  HECl-ADAPT,  this  Is  accomplished 
by  defining  each  stream  link  and  the  triangles  that  drain  to  It 
as  the  equivalent  HECl-ADAPT  sub-basin. 

The  triangles  that  make  up  each  sub -has In  are  identified 
through  a  computer  analysis  of  the  terrain  and  topologic 
information  stored  within  the  data  base.  Since  each  ADAPT  sub- 
basin  Is  typically  composed  of  more  than  one  triangle,  a  certain 
amount  of  lumping  (averaging)  must  occur  to  derive  the  parameters 
tdiich  characterize  the  sub-basin.  The  degree  of  lumping  Is  a 
function  of  triangle  size. 

ADAPT  also  has  the  capablll^  of  defining  sub -water sheds. 

Sub -watersheds  are  portions  of  the  overall  data  base  and  are 
identified  using  stream  links.  The  sub-watershed  definition  is  a 
windowing  capability  which  makes  It  possible  to  model  specific 
portions  of  the  data  base  without  having  to  use  the  entire  data 
base. 

2.5.2  Overland  Flow  Parameters 

Overland  flow  is  controlled  by  the  quantity  and  temporal 
pattern  of  rainfall,  by  infiltration  and  evaporative  losses,  and 
by  the  process  by  tdilch  water  travels  to  the  stream  channel. 

In  HECAD,  the  gage  locations  used  in  calculating  rainfall 


are  stored  in  the  ralngage  file.  Total  rainfall  is  computed  for 
each  sub-basin  using  a  weighting  function  in  which  the  weight  is 
inversely  proportional  to  the  distance  from  the  gage  to  the 
centroid  of  the  sub-basin.  The  user  can  limit  the  nuBd>er  of 
gages  used  for  each  sub -basin  by  specifying  a  maximum  distance  or 
a  maximum  nuiri>er  of  gages.  Up  to  five  gages  may  be  used  to 
define  storm  totals.  The  gage  closest  to  a  sub-basin  is  used  to 
define  the  temporal  distribution. 

Two  methods  are  available  for  modeling  losses  in  the  HECl- 
ADAPT  system:  initial  and  \miform  loss  rate  and  SCS  curve 
nuid>er  (CH)  (29).  Numerical  values  of  initial  and  uniform,  loss  rate 
and  CM  are  stored  within  the  calibration  file  as  a  function  of 
both  land  use  and  soil  hydrologic  group.  The  mix  of  land  use  and 
hydrologic  soil  group  within  each  triangle  is  stored  in  the  data 
base.  HECAD  first  calculates  loss  rate  parameters  on  a  triangle - 
by- triangle  basis  and  then  computes  an  areally -weighted  average 
for  each  sub-basin  by  accessing  the  Information  stored  in  both 
the  calibration  file  and  the  data  base. 

The  interface  also  Includes  an  option  for  adjusting  CN  based 
on  antecedent  precipitation  and  season. 

Both  methods  of  calculating  losses  apply  only  to  pervious 
areas.  The  Impervious  area  of  each  triangle  is  calculated  by 
HECAD  using  calibration  file  data  relating  land  use  to  percent 
imperviousness  and  land  use  data  from  the  data  base. 

Four  methods  are  available  in  HECl-ADAFT  for  transforming 
rainfall  excess  into  stib-basin  runoff  hydrographs:  Clark  Unit 


Graph,  Snyder  Unit  Graph,  SCS  Dlaenslonless  Graph  and  Klneoaclc 
Wave.  For  the  Clark  Unit  Graph  nethod,  HECAD  coapuces  a  Clne- 
area  curve  and  the  tuo  parameters:  time  of  concentration  and 
storage  coefficient.  Manning"  s  eqtiation  is  used  to  derive 
velocity  from  which  travel  time  is  computed.  The  roughness 
coefficient  is  derived  for  each  triangle  based  on  the  land  use  as 
stored  in  the  data  base  and  the  rou^ness  supplied  by  the 
calibration  file.  Slope  data  are  derived  from  the  data  base. 

The  travel  time  and  area  associated  with  each  triangular  element 
.'are  used  to-develop  a  time-area  curve  for  each  sub-basin.  Time 
of  concentration  and  storage  coefficient  are  computed 
respectively  as  the  longest  triangle  travel  time  for  the  sub- 
basin  and  as  the  areally-weigihted  travel  time  for  the  sub-basin. 
The  storage  coefficient  can  also  be  coaq>uted  by  specifying  a 
ratio  R/(Tc4fi.),  idiere  R  is  storage  coefficient  and  Tc  is  time  of 
concentration  for  use  with  all  sub-basins. 

Similar  procedures  are  used  to  derive  the  Snyder  and  SCS 
parameters.  For  the  Snyder  Method,  lag  is  computed  as  area- 
weighted  travel  time,  the  peaking  coefficient  is  supplied  as 
input  by  the  user  and  the  time-area  curve  is  developed  in  the 

same  way  as  for  the  Clark.  For  the  SCS  method,  lag  is  also 

$ 

calculated  as  the  area -weighted  travel  time. 

The  Kinematic  Wave  Method  parameters  are  slope,  roughness 
and  overland  flow  length.  Area-wei^ted  values  of  slope  and 
roughness  are  derived  using  the  slope  and  roug^ess  data  from  the 
data  base  and  roug^ess  data  stored  as  a  function  of  land  use  in 
the  calibration  file.  The  overland  flow  length  is  calculated  by 


one  Mthod  if  tuo  overland  flow  planes  converge  Co  a  cencral 
channel  and  by  a  different  neChod  when  a  single  overland  flow 
plane  drains  Co  Che  channel « 

2.5.3  SCrean  ParaaeCers 

The  three  stream  routing  nethods  included  in  HECl-ADAFT  are: 
Kineoatic  Wave,  Muskingua  and  noraal- depth.  All  Kinematic  Wave 
parameters  (channel  length,  slope,  roughness  shape,  width  and 
side  slopes)  are  calculated  or  extracted  from  the  drainage 
network  file  and  auxiliary  network  file.  Muskingum  K  is  asstimed 
to  be  equal  to  Che  reach  travel  time  'as  computed  using  Manning's 
equation.  The  n\ad>er  of  routing  steps  Is  computed  as  travel  time 
divided  by  the  time  step  parameter  supplied  as  input  by  the  user. 
Muskingum  X  is  also  supplied  by  the  \iser  and  Is  Che  same  for  all 
channels.  For  normal  depth  routing,  cross-section  data  and 
Manning's  roug^mess  are  stored  in  the  auxiliary  network  file. 
Reach  length,  slope  and  the  datum  elevation  are  derived  from  the 
data  base. 

Reservoir  routing,  base  flow  and  channel  loss  parameters  may 
also  be  input.  Reservoir  routing  is  accomplished  by  passing  flow 
throu^  links  identified  as  reservoirs  with  no  transformation. 
Storage  r<  ating  is  performed  only  when  the  downstream  end  of  a 
reservoir  link  is  encountered.  Parameters  required  for  reservoir 
routing  are  stored  in  the  reservoir  file. 

Base  flow  parameters  entered  by  the  user  are;  base  flow 
yield  in  cfs/sq\iare  mile,  a  base  flow  ratio  by  which  the  peak 
flow  is  multiplied  Co  determine  >dien  the  recession  part  of  the 
hydrograph  starts  and  a  recession  coefficient  describing  the 


slop*  of  the  rocesslon  curv*.  A  conscsnt  channel  loss  rate  (In 
cfs)  and  a  paraaeter  representing  the  percentage  of  reaalnlng 
flow  after  constant  loss  is  subtracted  out  nay  also  be  entered  by 
the  user. 

2.5.4  Plotting  Capability 

The  graphics  capability  of  ADAPT  is  one  of  the  systesis  most 
useful  aspects.  A  series  of  plots  used  to  develop  the  Castro 
Valley  nbdel  are  shorn  to  denonstrate  this  capability. 

Figure  1  shows  a  TIN  for  the  Castro  Valley.  This  plot  Is 
used  to  check  for  errors  In  the  vertex  and  triangle  files. 

Contour  and  slope  direction  plots.  Figures  2  and  3,  are  also  used 
for  error  checking.  The  contour  plot  can  be  overlayed  on  the 
original  topographic  map  to  spot-check  elevations.  The  slope 
direction  plot  is  useful  for  insuring  that  all  triangles  drain 
inward  along  the  watershed  boundaries.  If  the  topography  is  not 
modeled  satisfactorily,  the  triangle  network  may  require 
modification.  Figures  4  and  5  show  land  use  and  soil  polygons 
^ile  Figure  6  shows  the  unique  attribute  polygons.  The  last 
plot  of  this  group.  Figure  7,  shows  a  drainage  network  plot 
superimposed  on  the  contour  plot  and  dearanstrates  the  overlay 
capability  of  the  system.  Because  of  the  high  quality  oT  these 
plots,  any  of  them  can  be  used  as  figures  or  displays  In  reports. 
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3.  UTERATORE  REVIEW 


This  litsrature  review  discusses  papers  on  a  wide  range  of 
water  resource  applications  for  which  CIS's  have  been  developed. 
Papers  are  listed  in  chronological  order  for  ease  of  presentation 
and  to  provide  a  sense  of  the  technological  evolution  this  nethod 
has  experienced.  Sone  of  the  early  grid>based  systems  were  used 
strictly  as  organizational  aids.  The  new  TIN -based  systems  are 
being  employed  in  increasingly  sophisticated  models  which  take 
great  advantage  of  their  many  capabilites. 

Pentland  and  Cuthbert,  1971  (28).  This  paper  describes  a 
square  grid  nethod  used  to  automate  the  determination  of  regional 
hydrologic  relationships.  The  grid  method  provides  an  efficient 
means  of  integrating  hydrometric,  meteorologic  and  physiographic 
data.  Regression  analysis  is  used  to  define  mean  annual 
precipitation,  temperature  and  runoff  in  each  grid.  These 
results  and  physiographic  grid  data  are  used  in  a  second 
regression  analysis  to  define  monthly  flows  at  ungaged  sites.  A 
stochastic  model  is  then  applied  to  generate  synthetic  flows  for 
operational  hydrology. 

Huggins,  Burney,  Kunder  and  Monk,  1973  (22).  The  watershed 
model  described  in  this  paper  is  based  on  subdividing  catchments 
into  grids  idiich  are  assumed  uniform  with  respect  to  hydrologic 
variables.  Response  of  each  grid  is  characterized  by 
deterministic  equations.  Interaction  between  Individual  grid 
elements  and  conqrasite  watershed  response  is  analyzed  by 
integrating  the  continuity  of  mass  equation  over  the  ^diole  basin. 


SMdar.  1974  (32).  A  aodel  calUd  -OYIAM  II*  Is  used  Co 
project  land  use  pattems  for  Che  purpose  of  predicting  surface 
runoff.  Alcemacive  future  scenarios  are  investigated  Co  derive 
a  range  of  future  conditions.  A  grid  SMthod  is  used  to  iiqmt  and 
output  data. 

Grayoan,  Males,  Gates  and  Madder .  1975  (11).  This  paper 
describes  ADAPT  and  an  application  of  ADAPT  to  water  quality 
Modeling. 

Charbonneau.  Fortin  and  Morin,  1975  (1) .  The  Cequeau  -Model 
uses  a  grid  system  to  define  surface  elements.  The  model  assigns 
each  cell  a  maxiimm  and  minimum  teiqteraCure  and  precipitation 
value  based  on  computations  using  data  from  existing 
meteorological  stations.  A  hydrologic  balance  is  done  with 
"more  or  less*  sophisticated  math  models  which  describe 
individual  hydrologic  processes.  This  is  the  *production*  part 
of  the  model.  A  "transfer*  part  models  the  movement  of  water 
from  cell  to  cell  using  a  "transfer  coefficient*  which  is  a 
function  of  the  physiographic  characteristics  of  each  grid. 

HEC,  1975  (14).  This  report  describes  and  illustrates  the 

application  of  data  management  and  analytical  techniques 

» 

developed  by  the  HEC  for  application  in  comprehensive  flood  plain 
information  studies.  The  technique  uses  gridded  geographic  data 
to  analize  the  effects  of  alternate  land  use  pattems  on  flood 
hazard,  general  damage  potential  and  environmental  status  of  the 
study  area.  AUTOMAP  II,  a  program  developed  by  the  Environmental 
Systeau  Research  Institute  in  Redlands,  California,  is  used  to 


■wnlpulaCa  the  data  and  is  a  key  to  the  techniques  developed. 

Fabos  and  Joyner.  1976  (10).  The  nodel  in  this  paper 
provides  a  procedure  to  assess  special  resources  hazards  and 
developnent- suitability  potentials  to  aid  in  planning.  A 
mapping  system  called  "COMLUP"  is  used  to  develop  overlays  to 
form  composite  special  resource,  hazard  or  development 
suitability  maps.  The  mapping  system  utilizes  a  polygon  format 
for  inputting  data.  The  program  automatically  converts  from  this 
format  to  a  grid  representation  for  data  manipulation. 

Li,  Shahholtz,  Contractor  and  Carr.  1977  (25).  This  model 
involves  discretization  of  a  drainage  basin  into  hydrologic 
response  tmlts  (HRU's)  based  on  soils,  land  use,  and  physiographic 
features.  Precipitation  excess  is  generated  using  the  Holtan 
eqtiation  and  flow  is  routed  using  a  finite  element  solution  of 
the  kinematic  wave  equations.  A  grid  method  and  digital 
processing  are  used  to  derive  HRU's  from  overlays  of  soil  and 
land  use  naps.  Finite  elements  and  HRU  boundaries  do  not 
coincide  and  the  program  must  therefore  derive  a  weighted 
precipitaion  excess  for  each  element  based  on  HRU's  within  the 
element  before  flow  routing  is  done. 

Gupta  and  Solomon,  1977  (13).  In  this  model,  a  basin  is 
conceptualized  as  being  composed  of  a  set  of  finite-sized  grids 
with  each  grid  homogeneotis  in  physical  characters  tics.  The  data 
base  contains  a  series  of  digitized  maps  of  physiographic  data, 
time  series  data,  and  location  of  meteorological  and  hydrologic 
stations.  Kap  data  are  digitized  by  tislng  a  polygon  method.  A 
series  of  computer  programs  transpose  this  data  into  grid  data. 


Rainfall  excess  Is  coaputed  for  each  grid  using  Che  HolCan 
equacion.  A  surface  and  stib- surface  water  balance  is  done,  and 
flow  is  routed  using  the  Muskingum  oMthod. 

Davis,  1978  (2).  Spatial  data  management  techniques  for 
comprehensive  flood-plain  studies  ere  described  in  this  paper.  A 
grid  cell  format  is  used  to  score  data  such  as  existing  and 
future  land  use,  physiographic  data,  hydrologic  sub-basins,  and 
environmental  habitats.  UtiliQr  progrsms  access  files  and  create 
input  for  programs  used  In  flood  hazard  evaluation,  flood  damage 
analysis  and  environmental  assessments. 

Jett,  Weeks  and  Grayman.  1979  (23).  This  paper  describes  an 

application  of  Che  ADAPT  triangular  data  base  to  hydrologic 

modeling  using  several  alternative  rainfall -runoff  models.  The 

paper  emphasizes  that  since  GIS*s  provide  detailed  physical 

modeling  of  drainage  basins,  an  analysis  doesn’t  have  to  be 

constrained  to  acquiring  data  for  a  specific  model.  Instead,  one 

can  select  the  most  appropriate  model  based  on  Che  type  of 

investigation.  Hydrologic  models  developed  for  use  with  ADAPT 

range  from  simple  unit  hydrograph  models  with  unit  hydrograph 

parameters  derived  from  average  basin  characteristics,  to 
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detailed  routing  models  which  compute  excess  for  each  triangular 
element  and  route  flow  through  both  overland  and  stream  networks. 

Thomsen  and  Striffler,  1980  (33).  This  report  describes  a 
watershed  information  system  which  is  used  to  continuously 
simulate  snoirpack  processes  and  to  generate  stream  flow 
forecasts.  The  system  utilizes  remote  sensing  data  to 


periodically  update  the  simulatioa.  A  grid  approach  is  used  to 
create  a  set  of  overlays  containing  data  on  elevation,  aspect, 
vegetation  and  soils.  Two  prograas  derive  parameter  decks  for 
water  yield  and  stream  flow  oradels  using  these  overlays.  The  - 
water  yield  program  does  water  balance,  snow  accumulation  and 
melt  calculations.  Output  from  this  model  drives  the  stream  flow 
model  which  uses  a  Darcy* type  equation  and  the  continuity 
equation  to  calculate  lateral  flow.  Deep  seepage  and  baseflow 
are  treated  empirically.  The  model  does  not  consider  Hortonian- 
type  infiltration  because  infiltration  rates  on  terrain  simulated 
by  the  model  are  generally  much  greater  than  any  snowmelt  or  rain 
event. 

Eli,  Palmer  and  Hamrio,  1980  (5).  This  paper  describes  an 
application  of  ADAPT  to  high  resolution  modeling  of  an  abandoned 
strip  mine  in  Vest  Virginia.  The  model  consists  of  270  triangles 
some  of  which  are  a  fraction  of  an  acre.  The  object  of  the  study 
is  to  model  the  micro -topography  of  the  site  including  spoil 
piles,  access  roads,  benches  and  drainage  courses.  The  paper 
demonstrates  how  ADAPT  can  be  an  efficient  method  for  increasing 
hydrologic  model  resolution.  It  also  demonstrates  how  this 
increased  resolution  allows  accurate  modeling  of  flow  direction 
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and  concentration  of  runoff. 

Eli,  1981  (7).  This  paper  proposes  a  combination  of  ADAPT 
with  the  Hewlett  concept  of  variable  source  areas  of  runoff  for 
Continuous  or  single  event  modeling  on  small  watersheds.  The 
paper  describes  the  previous  application  of  ADAPT  for  surface 
mine  hydrology  in  West  Virginia  and  suggests  modifications  to 


original  roucines  that  will  enable  continuous  nodeling.  It  also 

outlines  Bodiflcation  of  the  overland  routing  schene  to 

incorporate  the  concept  of  contributing  area.  A  series  of 

"runoff  bands"  which  bound  the  contributing  area  are  detemined 

by  a  new  set  of  decision  rules.  Ihe  paper  describes  how  the 

runoff  bands  can  be  utilized  in  erosion  and  sediment  yield 

computations.  It  also  describes  how  below  ground  surfaces  can  be 

represented  by  assigning  smre  than  one  elevation  to  each  triangle 

vertex,  liiese  additional  surfaces  can  be  used  to  do  mass*  balance 

.  •  , 

for  contlt^us  hydrologic  modeling. 

Ell  and  Paulin,  1981  (6).  This  paper  describes  applications 
of  TIN  type  GIS*s  to  runoff  and  erosion* sedimentation  modeling. 

It  demonstrates  how  GIS's  can  be  used  to  derive  input  for 
existing  hydrologic  models  such  as  SCS  TR*20.  It  also  suggests 
an  alternative  method  for  computing  overland  flow  lengths. 

Instead  of  constructing  centroid- to-centroid  connecting  lines, 
the  downs  lope  vector  becomes  the  actual  flow  path.  Flow 
direction  changes  as  triangle  boundaries  are  crossed  and  triangle 
slopes  change.  The  paper  suggests  that  present  applications  do 
not  take  advantage  of  the  spatial  resolution  available  in  these 
models.  It  recommends  using  the  principles  developed  for 
"cascading  planes"  to  develop  a  flow  model  which  is  more 
cosq>atible  with  GIS's.  The  paper  also  mentions  a  microcomputer 
compatible  TIN  GIS  called  "GEOSPHERE"  which  is  being  developed  by 
Ell  for  small -watershed,  hlgh-resolutlon  environments. 

La  Garde,  1982  (24).  This  report  describes  a  rainfall- 


runoff  ■odol  idilch  onploys  «  polygon- Co*grld  CIS  Co  score  dai:s. 
The  reporc  provides  scep-by-scep  insCruccions  for  cresClng  a  daCa 
base  and  running  Che  nodal.  Soil,  land  use  and  Copographic  daCa 
are  inpuC  using  a  polygon  nethod.  Auxiliary  programs  converc  and 
process  Che  data  InCo  grid  fonaac  Co  creace  a  CIS.  Rainfall 
excess  Is  coiq>uCed  for  each  grid  using  CN's.  A  lag  equaclon 
which  Is  a  funcclon  of  CN  and  surface  slope  is  used  co  allocaCe 
flow  becween  grids.  Lag  divided  by  Cine  sCep  defines  Che 
fracClon  of  flow  In  Cenporary  sCorage  chac  will  be  removed  co  Che 
lowesc  downscreaa  grid.  A  flow  hisCory  can  be  developed  for  each 
grid. 

Graynan,  Kales,  GaCes  and  Harris,  1982  (12).  This  paper 
describes  appllcaClons  of  AJDAPT  Co  urban  hydrology.  The 
advancages  of  ADAN  for  modeling  urban  hydrology  are  oucllned 
Including  its  ability  to  provide  a  continuous  model  of 
topography.  The  paper  describes  how  ADAPT  can  be  used  Co  model 
both  natural  and  man-made  networks.  It  illustrates  application 
of  ADAPT  Co  urban  hydrology  with  example  projects  in  Ohio, 

Vtyomlng  and  Pennsylvania.  The  Ohio  study  involved  detailed 
ralnfall-runoff/non-polnc  source  pollution  modeling  of  twelve 
nortlieasCem  Ohio  sub-basins.  In  the  Vyomlng  study,  the,  issue 
was  determination  of  the  Impact  of  proposed  future  development  on 
an  existing  sewer  and  drainage  system.  In  Pennsylvania,  the 
stuity  involved  rainfall-runoff  modeling  and  generation  of  flood 
plain  maps  for  Che  main  stream  drainages. 

Ell  and  Paulin,  1983  (8).  This  paper  describes  a 
sensitivity  analysis  of  a  rainfall -ninoff  model  consisting  of  the 


ADAPT  syscen  and  a  linear  res-.  .r>llnear  channel  routing  «odel. 
CN  Is  used  to  generate  excess  precipitation.  The  nuaber  of 
triangular  eleaents  is  varied  to  test  sensitivity  of  the  outflow 
hydrograph  to  terrain  sK>del  resolution.  Lag-coefficients  are 
also  modified  to  determine  the  effect  on  model  results.  Using 
three  different  model  resolutions  of  the  basin,  it  is 
demonstrated  that  as  the  number  of  triangles  is  increased, 
average  link  slopes  Increase,  maximum  triangle  slopes  increase, 
average  triangle  areas  decrease,  and  number  of  stream  links 
increase. 

It  is  also  demonstrated  that  model  results  are  a  function  of 
lag  coefficients  chosen.  In  cases  where  the  proportion  of  lag 
assigned  to  the  linear  reservoir  is  50  percent  or  less,  the  high 
resolution  model  peaks  sooner  and  higher  than  the  low  resolution 
model.  The  situation  reverses  when  more  than  50  percent  of  the 
lag  is  assigned  to  the  linear  channel.  It  is  concluded  that  for 
"realistic"  values  of  the  lag  coefficients,  the  model  does  not 
require  a  high  resolution  representation  to  yield  acceptable 
results . 

Heggen,  1983  (20).  This  CIS  employs  a  grid  representation 
of  the  watershed.  Each  grid  is  described  by  elevation,’  soil  and 
cover  characteristics,  and  channel  descriptions  if  applicable. 

The  CN  method  is  used  to  define  surface  infiltration.  Channel 
infiltration  (important  in  New  Mexico)  is  estimated  using  an 
empirical  expression  developed  for  New  Mexico.  Surface  nmoff  is 
described  using  Manning's  equation.  Effective  slope  and  slope 
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l«ngth  ar*  darivad  by  aaplrical  ralacionshlpa  and  by  flald 
astlaatas  raapactlvaly.  Iha  dlrecclon  of  channal  flow  is 
coaputad  by  a  partitioning  routine  which  divides  outflow  by  grid 
based  on  relative  grid  elevations.  A  set  of  channel  hydraulic, 
characteristics  nnist  be  assuaed  to  accomplish  this. 

Eli,  1983  (9).  This  paper  describes  the  application  of 
GIS's  to  planning,  design  and  analysis  of  coal  mines.  It 
presents  an  overview  of  available  GIS's  and  discusses  advantages 
and  disadvantages  of  each.  A  description  of  a  new  TIN  based 
system  called  "HYGIS"  (Hybrid  Geographic  Information  System)  is 
presented  in  the  paper.  This  system  uses  a  TIN  to  represent 
three-dimensional  surfaces  above  and  below  ground  Two- 
dimensional  polygon  overlays  containing  attribute  information  can 
be  created  independent  of  the  TIN'S.  A  three-dimensional  grid 
cell  system  is  incorporated  to  aid  in  locating  specific  areas  of 
the  data  base.  Grid  cell  structure  also  aids  in  connecting  the 
multiple  TIN  surfaces  and  overlays.  The  system  is  used  to 
produce  various  maps,  including  projections  and  cross-sections  of 
surface  and  s\ib- surface  structures.  Engineering  data  including 
lengths,  areas  and  volumes  can  also  be  calculated. 

HEC,  1983  (16).  This  document  describes  the  procedure  for 
developing  HEC-1  input  data  using  a  grid  cell  GIS.  Data  is 
entered  in  a  grid  format  using  a  program  called  "BANK." 
Verification  of  input  is  accomplished  with  program  "RIA"  which 
displays  stored  data  using  line  printer  graphics.  Program 
"HYDPAR"  is  the  interface  between  the  grid  cell  GIS  and  the 
rainfall -runoff  model  (HEC-1) .  HYDPAR  derives  loss  rate  and  unit 


hydrograph  paranecers  fron  the  CIS.  Results  are  output  to  a  file 
which  can  be  autonatlcally  transferred  to  HEC-1.  SCS  CN  and 
percent  iatpervlousness  are  derived  using  HYDPAR  as  are  the  SCS 
and  Snyder  unit  graph  coefficients.  The  SCS  unit  graph  lag  is 
cooputed  using  an  equation  in  %rhich  the  lag  is  a  function  of 
average  basin  slope  and  CN.  Slope  and  CN  are  input  for  each 
grid.  Snyder's  lag  is  a  function  of  stream  lengths,  stream  slope 
and  percent  Imperviousness.  All  these  values  must  be  manually 
derived  and  input  to  run  HYDPAR. 

McKim,  Unger,  Merry  and  Ganthier.  1984  (26).  The  objective 
of  this  study  was  to  integrate  remotely  sensed  land  cover  data 
with  a  hydrologic  model  developed  for  the  Saginaw  River  Basin  in 
Michigan.  The  data  base  developed  was  compatible  with  the  HEC 
Spatial  Analysis  Methodology  (HEC-SAM)  software  (2) .  Two 
computer  programs  were  used  to  classify  land  use  from  the  Landsat 
images.  The  resulting  1.1  acre  Landsat  land  cover  classification 
was  converted  to  40>acre  grid  cells  using  an  aggregation  scheme. 
HEC-1  optimization  methods  were  used  to  derive  Clark  and  Snyder 
\mit  graph  parameters.  For  seven  gaged  sub -basins,  multiple 
linear  regression  was  then  used  to  develop  relationships  between 
unit  graph  parameters  and  the  land  use  classification  for  each 
sub-basin. 

Hong  and  Ell,  1985  (21).  This  paper  describes  a  rainfall- 
runoff  model  which  accounts  for  both  the  overland  flow- interflow 
and  the  infiltration- exfiltration  processes.  The  model  uses  a 
TIM- type  topographical  model.  Flow  direction,  slope,  hydrologic 
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and  topographic  characteristics  are  stored  In  the  DIM.  Using 
this  Infor.utlon,  a  progran  detemlnes  Che  series  of  elements 
which  contribute  to  each  stream  segment.  Each  series  is  treated 
as  a  set  of  planes  over  or  through  which  flow  passes.  Water  is 
routed  continuously  through  a  cos^ination  of  overland  flow  and 
interflow  from  the  top  element  down  to  the  stream  segment. 
Kinematic  Wave  routing  is  iised  to  describe  overland  flow  routing 
%ihile  Darcy's  law  is  used  to  describe  interflow.  The  storage* 
discharge  history  of  each  element  is  based  on  conservation  of 


mass. 


4.  TESTll^  PROGRAM 


This  secClon  of  the  report  describes  the  procedures  and 
results  of  the  testing  progran.  As  nentloned  previously  in  the. 
objectives  section,  the  aajor  goals  of  this  study  are:  1)  to  test 
the  ability  of  HECl -ADAPT  to  model  rainfall -runoff  on  ungaged 
basins;  2)  to  determine  the  effect  of  model  resolution  on  the 
simulated  outflow  hy  rograph;  and  3)  to  test  the  flexibility  of 
the  HECl-ADAPT  system. 

To  acc<»pllsh  these  goals,  two  drainage  basins  are  used  In 
the  testing  program.  The  first  is  Castro  Valley,  a  predominately 
urban  basin  of  5.5  square  miles  located  In  the  San  Francisco  Bay 
area.  Potter  Valley,  the  second,  is  an  agricultural  basin  with 
an  area  of  92.2  square  miles  located  in  the  Russian  River  basin 
In  northern  California.  These  basins  were  chosen  because  they 
represent  a  fairly  wide  range  of  geographic  conditions.  Modeling 
these  two  basins  should  provide  a  good  test  of  the  flexibility 
and  robustness  of  the  HECl-ADAPT  methodology  thus  accomplishing 
the  third  goal  of  the  study.  Castro  Valley  Is  modeled  using  two 
resolutions  to  accomplish  the  second  goal  of  the  study.  Potter 
Valley  Is  modeled  using  one  resolution. 

Both  basins  are  first  modeled  as  If  they  are  ungaged  to 
accomplish  the  first  goal  of  the  stu<ty.  Results  of  this  modeling 
are  highly  dependent  on  the  adopted  model  parameters.  Thus,  it  is 
iiq>ortant  to  chose  the  appropriate  curve  nun^er  to  use  with  a 
given  combination  of  land  use  and  soil  ^pe  and  the  appropriate 
roughness  and  percent  iiiq>ervlousness  to  associate  with  a  given 
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land  use.  Results  of  Che  unpaged  aodeling  effort  are  compared 
with  observed  data. 

The  higher  resolution  Castro  Valley  nodel  and  the  Potter 
Valley  sKMSel  are  then  calibrated  using  several  observed  flood 
events.  Results  of  the  calibrations  are  compared  with  historical 
data  and  with  hydrographs  generated  using  Clark  unit  graphs 
derived  by  HEC-1  optimization  methods.  The  models  are  Chen 
validated  using  other  historical  flood  events.  A  sensitivity 
analysis  is  performed  on  the  calibrated  models  of  the  two-  basins . 
Results  are  tabulated  and  discussed. 

Lastly,  modifications  to  Che  models  are  examined  and  some 
preliminary  runs  are  used  to  demonstrate  the  effect  of  these 
changes  on  simulation  results.  These  results  are  then  discussed. 


S.  AHALYSIS  AND  RESULTS 


5.1  Castro  Valley 

The  representation  of  topography  in  the  HECl -ADAPT  system  is 
probably  the  most  important  feature  of  the  method  since  derivation 
of  all  the  HEC-1  input  is  dependent  on  it.  In  this  testing 
program,  the  DIM  is  not  modified  to  incorporate  the  land  use  and 
soil  type  boundaries.  Instead,  this  information  is  input  on  a 
triangle -by- triangle  basis  with  each  triangle  containing  a  mix  of 
the  various  soils  and  land  uses  (see  Figure  8).  The  interface 
program,  HECAD,  determines  a  weighted  average  land  use  and  soil 
type  for  each  triangle  based  on  the  percentages  in  each  triangle. 
This  approach  is  used  throu^out  the  testing  program. 

Castro  Valley  was  the  first  basin  modeled  using  HECl -ADAPT. 
This  basin  was  used  because  of  its  small  size  (5.5  square  miles) 
and  because  of  the  availability  of  all  necessary  data  at  the 
offices  of  the  Hydrologic  Engineering  Center.  Two  models  of 
Castro  Valley  were  developed.  The  first  model  consisted  of  39 
triangles  (see  Figure  9).  Average  triangle  area  for  this  model 
was  about  90  acres. 

Because  each  of  the  laodels  in  this  testing  program  is  first 
developed  assuming  the  basins  are  ungaged,  the  appropriate 
methods  for  computing  and  transforming  rainfall  excess  are 
dependent  on  available  data  and  on  basin  characteristics.  HECl- 
ADAPT  provides  two  methods  for  computing  rainfall  excess:  the 
SCS  QT  method,  and  the  inltial/uxtlform  loss  method.  In  the 
uncalibrated  models  developed  for  this  testing  program,  CN's  are 
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used  to  aodol  excess  because  CH's  can  be  related  to  soil  type  and 
land  use.  A  CN  adjustaent  subroutine  autoaatically  adjusts  CN's 
based  on  the  season  and  antecedent  precipitation  entered  by  the 
user.  Since  the  Castro  Valley  is  predominately  urban,  the 
Kinematic  Vave  model  is  used  for  both  overland  and  channel 
routing.  The  land  use  and  soil  data  tised  to  model  Castro  Valley 
were  taken  from  a  previous  HEC  study  of  Castro  Valley  (19). 

Table  1  tabulates  land  use,  soil  type,  percent  imperviousness  and 
CN  tised  in  the  uncalibrated  model  of  Castro  Valley. 

Given  the  geographic  information  stored  within  the  model  and 
the  calibration  data  input  by  the  user,  HECAD  derives  the  HEC-1 
omdel  coefficients  and  generates  HEC-1  input  data.  The  HEC-1 
input  data  generated  by  HECAD  for  the  39 -triangle  Castro  Valley 
model  are  on  file  at  the  HEC. 

HEC-1  was  run  using  this  input  data.  Figure  10  shows  the 
computed  and  observed  hydrographs  for  the  Jan  16,  1973  storm 
event.  It  can  be  seen  by  comparing  these  two  hydrographs  that 
the  observed  hydrograph  peaks  sooner  and  is  quite  a  bit  more 
peaked  than  the  computed  hydrograph.  The  observed  hydrograph 
peak  is  about  2.3  hours  before  the  computed  and  is  about  13 
percent  greater.  Rimoff  volumes  are  similar.  Table  2  gives  a 
tabulated  comparison  of  computed  and  observed  hydrographs. 

The  greatest  difference  is  in  the  hydrograph  timing.  Many 
things  could  be  affecting  the  timing.  For  example,  the  model  may 
not  adequately  represent  the  basin,  the  input  parameters  may  be 
inappropriate,  or  the  temporal  and  areal  distribution  of 
precipitation  may  not  be  representative.  To  test  the  adequacy  of 


TABLE  1 


CASTRO  VALLEY  LAUD  USE  AMD  SOILS  DATA 


land  use  CURVE  NUMBER  PERCENT 

Hyi«OIOGIC  SOIL  CROUP  IMPERVIOUSNESS 


A 

B 

C 

D 

Natural  Vegetation 

39 

61 

74 

80 

0 

Low  Density 

57 

72 

81 

86 

30 

Residential 

Mediua  Density 

61 

75 

83 

87 

40 

Residential 

Hlgih  Density 

82 

88 

92 

94 

75 

R«s  ldential/C<Nnercial 


Source:  HEC,  tindeted.  Oconee  Style  Hydrology  Workshop.  Urban 
Hydrology  Cource  Workshop  for  Castro  Valley. 
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FIGURE  10  CASTRO  VALLEY  39-TRIANGLE  MODEL 

UNCALIBRATED  JANUARY  16,  1973  EVENT 


TABLE  2 


GOKPARISCM  OF  GOHPUTED  AND  OBSERVED  HYMOCTAPHS 
CASTRO  VALLEY  39-TRIAHGLE  MODEL 
UNCAUBRATED  JANUARY  16,  1973  EVENT 


SUM  OF 

BQUIV. 

MEAN 

TIME  TO 

PEAK 

TIME  TO 

FLOWS 

DEPTH 

FLOW 

CENTER 

OF  MASS 

FLOW 

PEAK 

(cfs-lOain) 

(in) 

(cfs) 

(hrs) 

(cfs) 

(hrs) 

ConpuCed 

Hydrograph 

16105 

0.749 

248 

6.22 

467 

6.67 

Observed 

Ifydrograph 

17877 

0.832 

275 

4.91 

537 

4.33 

DIFFERENCE 

-1772 

-0.082 

-27 

1.31 

-70 

2.33 

PERCENT 

DIFFERENCE 

-9.91 

-9.91 

-9.91 

26.72 

-13.08 

53.81 

Um  39-trlaiigle  aodel  and  to  invoatlsata  die  laqpact  of  using  a 
higher  resolution  nodel,  the  82 -triangle  Castro  Valley  iiodel  was 
developed.  This  second  Castro  Valley  atodel  is  described  in  the 
next  section. 

5.1.1  Castro  Valiev  82 -triangle  Model 

Bl^ty-two  triangles  were  used  In  the  second  Castro  Valley 
■odel  to  attain  a  higher  degree  of  accuracy  In  the  topographic 
representation.  Average  triangle  area  Is  about  45  acres  which  is 
about  half  that  of  the  39 -triangle  nodel.  Figure  11  shows  the 
82 -triangle  representation.  Figures' 12  and  13,  respectively, 
show  contour  plots  developed  for  the  39-  and  82- triangle  nodels 
using  the  graphics  capability  of  HECl-ADAPT.  One  can  see  by 
cooparlng  these  plots  that  the  82-trlangle  nodel  has  sone  steeper 
slopes.  This  nodel  also  adds  two  of  the  snaller  tributaries  to 
the  representation  of  the  channel  systen. 

Again,  the  basin  is  first  modeled  as  if  it  were  ungaged. 
Thus,  the  only  change  between  the  first  and  second  Castro  Valley 
models  is  the  topographic  representation.  The  calibration 
parameters  (CN's,  rou^nesses)  renain  the  sane  for  the  initial 
runs.  The  82- triangle  nodel  Is  later  calibrated  using  several 
historical  events.  The  39 -triangle  nodel  Is  not  calibrated. 

The  effect  of  this  higher  resolution  is  apparent  in  a 
conqiarison  of  the  two  hydrographs  computed  using  the  different 
nodels.  Figure  14  shows  the  39-  and  82-trlangle  model 
hydrographs  and  the  observed  hydrographs  for  the  January  16,  1973 
atom.  Table  3  gives  a  tabulated  coiq^arlson  of  the  two  nodels. 
The  82-trlangle  nodel  appears  to  concentrate  runoff  faster  than 
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TABLE  3 


COMPARISON  OF  39-  AND  82-TRlANGLE  CASTRO  VALLEY 
MODEL  HYDROGRAPHS  UHCALIBRATED 
JANUARY  16.  1973  EVENT 


SUM  OF 

BQUIV. 

MEAN 

TIME  TO 

PEAK 

TIME  TO 

FLOWS 

DEPTH 

FLOW 

CENTER 

OF  MASS 

FLOW 

PEAK 

(cfS'lOmln) 

(in) 

(cfs) 

(hrs) 

(cfs) 

(hrs) 

39 -Triangle 
Model 

16105 

0.749 

248 

6.22 

467 

6.67 

82-Trlangle 

Model 

16871 

0.763 

260 

6.03 

506 

4.67 

DIFFERENCE 

-766 

-0.014 

CM 

« 

0.19 

-39 

2.00 

PERCENT 

-4.54 

-1.83* 

-4.54 

3.15 

-7.71 

42.82 

DIFFERENCE 


*  Discrepancy  between  sua  of  flows  and  equivalent  depth  Is  caused 
by  slight  differences  In  drainage  area  between  the  .two  models. 


tile  39 -triangle  aodel  because  of  the  steeper  slopes  and  Increased 
nuiriier  of  stream  links.  The  82 -triangle  SKidel  hydrograph  has  a 
steeper  rising  limb  and  a  greater  peak  than  the  39 -triangle  model 
hydrograph.  The  lag  (center  of  auiss  to  center  of  mass)  of  the  - 
82-trlangle  model  is  also  less  than  that  for  the  39-triangle 
model.  Although  the  differences  are  not  great,  the  results 
appear  to  agree  with  those  obtained  by  Eli  (8) . 

The  82 -triangle  model  response  is. still  quite  a  bit  slower 
than  that  of  the  observed  basin  however.  This  is  evident  in  the 
slower  lag  time  and  the  smoothness  of  the  computed  hydrograph. 

The  observed  Castro  Valley  response  is  almost  immediate  as  shown 
by  the  rapidly  rising  and  falling  limbs  of  the  observed 
hydrograph.  There  is  apparently  little  basin  storage.  Because 
the  basin  is  predominately  urban,  much  of  it  is  drained  by 
gutters  and  storm  sewers  and  the  present  aodel  has  no  direct 
ability  to  account  for  this.  Since  tiie  volume  of  the  runoff  in 
the  simulation  run  is  comparable  to  the  observed,  calibration  of 
the  aodel  is  first  approached  by  lowering  roughness  factors  to 
get  a  quicker  response  to  compensate  for  the  unmodeled  storm 
drainage  system. 

The  iaqiact  of  modifying  channel  and  overland  roughnesses  is 
illustrated  on  Figure  IS  for  the  January  16,  1973  flood.  This 
model  responds  more  quickly  as  is  evident  in  the  steeper  rising 
limb,  the  smaller  lag,  and  the  spikiness  of  the  hydrograph.  One 
BKire  flood  event  is  simulated  using  this  aodel  and  the  results 
are  shown  on  Flgi^nre  16. 

The  ruiKiff  volumes  are  low  and  the  response  slow  in  the 


siaulaClons  of  diose  evencs.  To  obcoin  a  beCCer  callbraCion,  the 
initial/unlfoni  loss  rata  function  is  used  to  get  both  the 
correct  runoff  voluaws  and  to  shift  aore  of  the  runoff  voluae  to 
the  rising  liab  of  the  hydrograph.  Roughnesses  are  again 
adjusted  until  a  reasonable  natch  between  the  actual  and  coaqputed 
hydrograph  is  obtained.  The  HEC-1  input  data  for  the  calibrated 
nodel  are  on  file  at  the  HEC.  Results  of  this  calibration  are 
shown  in  Figures  17a  and  18a.  For  coiq>arison,  hydrographs 
generated  by  the  previous  oMdified  roughness  nodel  are  also 
shown.  Figures  17b  and  18b  coaq>are  tbe  calibrated  nodel 
hydrographs  Mixh  hydrographs  generated  by  HEC*1  using  optinlzed 
Clark  unit  Iqfdrograpli  paraneters  (19).  Table  4  gives  a  tabulated 
ecnqparlson  of  the  observed  and  calibrated  BK>del  hydrographs  for 
the  January  16,  1973  flood  event. 

For  the  calibrated  nodel,  runoff  volune  is  about  seven 
percent  less  than  the  observed  voluae  while  the  nodel  lag  (center 
of  nass  to  center  of  nass)  is  about  six  percent  greater  than  the 
observed  lag.  Model  peak  flow  was  about  28  percent  greater  than 
the  observed  peak  flow. 

One  can  see  from  Figures  17b  and  18b  that  the  HECl -ADAPT 
sinulations  ax^  the  HEC-1  sinulations  using  optinlzed  Clark 
paraneters  produce  hydrograph  peaks,  voluaes  and  tlning  that  are 
quite  similar.  After  conpletlng  the  calibration  runs,  both 
aodels  are  verified  using  the  Decefld>er  22,  1971  flood  event. 
Figure  19  shows  the  hydrographs  for  these  sinulations  while  Table 
S  conpares  the  simulation  results. 

AiriMMigh  both  models  do  a  poor  Job  of  reproducing  this 
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82.TRI  FIGURE  18a  CASTRO  VALLEY  82-TRIANGLE  MODEL 
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TABLE  4 


COMPARISON  OF  COMPUTED  AND  OBSERVED  HYlXtOGRAPHS 
CASTRO  VALLEY  82-TRIANGLE  MODEL 
CAUBBATED  JANUARY  16,  1973  EVENT 


SUM  OF 

EQUIV. 

MEAN 

TIME  TO 

PEAK 

FLOWS 

DEPTH 

FLOW 

CENTER 

OF  MASS 

FLOW 

(cfs-lOnln) 

(in) 

(cfs) 

(hrs) 

(cfs) 

CooqpuCed 

Hydrograph 

16648 

0.753 

256 

5.19 

686 

Obsarved 

Hydrograph 

17877 

0.809 

275 

4.91 

537 

DIFFERENCE 

-1229 

•0.056 

-19 

0.28 

149 

PERCENT 

-6.88 

-6.88 

•6.88 

5.69 

27.66 

DIFFERENCE 


TIME  TO 
PEAK 

(hrs) 

3.67 

4. .33.  • 

-0.67 

-15.47 
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TABLE  5 

OOMPARISON  OF  82-TRlAMGLB  MODEL  AND  HEC-1 
OPTIMIZED  CLABK  MODEL  HYDROGRAPHS 
DECEMBER  22.  1971  EVENT 


SUM  OF 

BCHIIV. 

MEAN 

TIME  TO 

PEAK 

TIME  TO 

FLOWS 

DEPTH 

FLOW 

CENTER 

OF  HASS 

FLOW 

PEAK 

(cfs'lOmln) 

(in) 

(cfs) 

(hrs) 

(cfs) 

(hrs) 

82 -Triangle 
Model 

2410 

0.109 

69 

3.85 

197 

3.00 

Optimized 
Clark  Model* 

2711 

0.143 

77 

3.76 

220 

2.83 

DIFFERENCE 

-301 

-0.034 

-8 

0.09 

-23 

0.17 

PERCENT 

DIFFERENCE 

-11.10 

-23.78+ 

-11.10  2.39 

-10.45 

6.01 

OBSERVED 

3204 

0.169 

92 

3.06 

580 

3.00 

*  Source:  HEC,  undeCed.  Oconee  SQrle  Hydrology  Workshop.  Urban 
Hydrology  Course  Workshop  for  Castro  Valley. 

-i-  Discrepancy  between  sum  of  flows  and  equivalent  depth  is  caused  by 
slight  differences  In  drainage  area  between  the  two  models. 


event,  possibly  due  to  innaccurate  strean  or  rain  gage  data,  the 
hydrographs  produced  by  the  two  itodels  are  quite  sinilar.  From 
Table  5  it  can  be  seen  that  simulated  hydrograph  parameters  only 
vary  by  about  10  percent.  HECl -ADAPT  appears  to  provide  a 
physically-based  methodology  which  simulates  the  rainfall -nmoff 
process  about  as  successfully  as  the  HEC-1  optimized  Clark 
method. 

5.1.2  Sensitivity  Analysis 

In  a  sensitivity  analysis,  the  major  point  that  needs  to  be 
established  is  the  relative  sensitivity  of  state  variables  (like 
peak  discharge)  to  changes  in  the  values  of  model  parameters 
(like  surface  roughness  factors).  The  results  of  a  sensitivity 
analysis  give  the  modeler  a  *feel"  for  the  effect  inaccurate 
parameter  estimation  may  have  on  the  simulation.  The  definition 
below  allows  the  iMdeler  to  decide  what  parameters  have  the  most 
and  the  least  impact  on  model  results  (27) . 

Sij-(ACi/Ci)/(ABj/Bj) 

Where:  Sij-sensitivity  coefficient 

ACi«^hange  in  state  variable  (e.g.  discharge) 

Ci  -reference  value  of  state  variable 
ABj-change  in  parameter  (e.g.  channel  roughness) 

Bj  -reference  value  of  parameter 

All  parameters  are  held  constant  except  the  one  being 
studied  to  isolate  the  impact  of  the  Individual  parameter.  Using 
this  approach,  the  sensitivity  coefficients  are  computed  and  then 
compared  directly  to  determine  what  parameters  have  the  greatest 
Impact  on  model  results.  Table  6  shows  the  results  of  the 
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TABLE  6 

SENSITIVITY  ANALYSIS 
CASTRO  VALLEY  82-TRIANGLE  MODEL 


PARAMETER 

CASE 

PEAK 

DISCH 

LAG 

EQUIV 

DEPTH 

SENSITIVITY  COEFFICIENTS 

aQ/q*  al/l* 

A  B/B  A  B/fi  A  B/B 

(%) 

(cfs) 

(hrs) 

(in) 

Overland 

Ref. 

686 

5.19 

0-753 

Roughness 

•1-20 

637 

5.33 

0.738 

0.357 

0.135 

0.100 

■HO 

661 

5.26 

0.745 

0.364 

0.135 

0.106 

-10 

717  * 

5.10 

0.760 

0.452 

0.173  • 

0.093 

-20 

746 

5.02 

0.768 

0.437 

0.164 

0.100 

Channel 

Ref. 

686 

5.19 

0.753 

Rouf^ess 

•t-20 

678 

5.20 

0.753 

0.058 

0.010 

0.0 

•flO 

682 

5.19 

0.753 

0.058 

0.0 

0.0 

-10 

690 

5.18 

0.753 

0.058 

0.019 

0.0 

-20 

694 

5.17 

0.753 

0.058 

0.019 

0.0 

Initial/ 

Ref. 

686 

5.19 

0.753 

Uniform 

•►20 

624 

5.23 

0.688 

0.452 

0.039 

0.432 

Loss 

•HO 

653 

5.21 

0.719 

0.481 

0.039 

0.452 

-10 

720 

5.16 

0.789 

0.496 

0.058 

0.478 

-20 

756 

5.14 

0.831 

0.510 

0.048 

0.518 

Percent 

Ref. 

686 

5.19 

0.753 

Impervious 

-►20 

729 

5.14 

0.812 

0.313 

0.048 

0.392 

HO 

705 

5.17 

0.778 

0.277 

0.039 

0.332 

-10 

667 

5.21 

0.724 

0.277 

0.039 

0.385 

-20 

645 

5.23 

0.695 

0.299 

0.039 

0.385 

*  AQ  »  change  in  peak  discharge 

Q  —  reference  value  of  peak  discharge 
AL  -  change  in  lag 
L  -  reference  value  of  lag 
AD  ■■  change  in  equivalent  depth 
D  “  reference  value  of  equivalent  depth 
AB  -  ci^ge  in  parameter  (e.g.  overland  rou^ness) 

B  «  reference  value  of  parameter 
AQA)  -  sensitivity  coefficient  for  peak  discharge 
AB/B 

Note;  reference  values  for  state  variables  and  parameters 
are  from  the  calibrated  model 


senslclvity  «naly»ls  for  the  82* triangle  Cestro  Valley  model. 

Prom  Table  6,  it  Is  seen  that  peak  discharge  and  equivalent  depth 
are  most  sensitive  to  estimates  of  loss  rates  while  lag  is  most 
sensitive  to  estimates  of  overland  roughness. 

Castro  Valley  responds  rapidly  to  rainfall  because  of  its 
comparatively  small  size  (5.S  square  miles)  and  its  urban 
character.  All  excess  produced  by  the  basin  is  transported  to 
the  basin  outlet  over  a  very  short  period  of  time.  This  rapid 
concentration  of  runoff  is  probably  the  reason  simulated  peak 
discharge  is  isost  sensitive  to  the  estimates  of  the  loss 
parameter. 

Peak  discharge  is  also  very  sensitive  to  estimates  of 
overland  roughness  because  of  the  direct  impact  of  this  parameter 
on  the  timing  of  runoff.  The  estiaute  of  channel  roughnesses  is 
less  significant  for  the  peak  for  two  reasons.  First,  the 
contribution  of  stream  travel  time  to  total  travel  time  is 
proportionately  less  for  small  basins.  Second,  the  channel 
roughness  factors  are  small  in  magnitude  to  begin  with.  A  10  or 
20  percent  change  in  a  suller  magnitude  parameter  will  not 
affect  the  simulation  as  much  as  a  10  or  20  percent  change  in  a 
larger  magnitude  parameter. 

Percent  imperviousness  has  a  smaller  impact  on  the  peaks 
becatise  the  percentage  of  lmpervio\is  surfaces  is  much  less  than 
the  percentage  of  pervious  stirfaces  for  the  basin. 

Lag  is  most  sensitive  to  estimates  of  the  overland  roughness 
parameters  because  of  the  direct  intact  of  this  parameter  on 
runoff  timing.  It  is  less  sensitive  to  the  estimates  of  channel 
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rou|^MSS  for  th«  saae  cwo  reasons  discussed  for  the  peak 
discharge. 

Equivalent  depth  is  nost  sensitive  to  the  estimates  of  loss 
parameters  and  percent  imperviousness  because  these  parameters 
determine  the  runoff  volume.  Equivalent  depth  is  most  sensitive 
to  the  loss  function  parameter  because  the  pervious  basin  area  is 
much  greater  than  the  impervious  basin  area.  Thus,  ninoff  from 
the  pervious  area  of  the  Castro  Valley  model  will  be  greater  than 
runoff  from  the  impervious  area. 

5.2  Potter  Valley 

Potter  Valley  was  the  second  basin  modeled  using  the  HECl* 
ADAPT  system.  This  basin  of  92.2  square  miles  is  much  larger 
than  Castro  Valley  and  is  mostly  woodlands  with  some  grasslands, 
cultivated  orchards  and  vineyards.  The  model  developed  for 
Potter  Valley  consisted  of  299  triangles.  Average  triangle  area 
was  about  200  acres  (see  Figure  20).  Soil  and  land  use  data 
were  obtained  from  the  U.S.  Soil  Conservation  Service  and  from 
the  California  Department  of  Water  Resources  respectively  (31,4). 
Table  7  tabulates  land  use,  soil  type  and  CN  for  Potter  Valley. 

The  modeling  of  Potter  Valley  serves  to  illustrate  a  major 
problem  that  must  be  addressed  when  using  HECl -ADAPT  on  larger 
basins.  The  topographic  model  of  Potter  Valley  is  good  where 
existing  channels  are  represented;  however,  the  topographic  model 
is  not  as  good  where  channels  exist,  but  are  not  modeled. 

Overlax>d  slopes  in  these  unmodeled  channel  areas  can  be  much  less 
than  the  actual  slopes.  The  dilemma  facing  the  user  is  deciding 


TABLE  7 


POTTER  VALLEY  LAHD  USE  AND  SOILS  DATA* 


I^AMD  use  CURVE  HUKBER* 

HYMtOLOGIC  SOIL  GROUP 


A++  B 

C 

D 

Chaparral 

44 

60 

66 

Grass -Oak 

46 

62 

67 

Irrigated  Pasture 

49 

65 

70 

Orchard 

53 

67 

71 

Woods -Forest 

55 

70 

77 

*  Source:  OUR,  X972.  Mendocino  CounQr  Land  Use  Maps. 

SCS,  1984.  Unpublished  Soil  Survey  Data  for 
eastern  Mendocino  County. 

+  Source:  SCS,  1972.  National  Engineering  Handbook,  Section  4, 
Hydrology. 

SCS,  1975.  Urban  Hydrology  for  Snail  Watersheds, 
Tecnical  Release  No.  55. 

++  No  "A*  soils  are  found  in  Potter  Valley. 


liiult  resolution  Is  necessary  to  give  reasonable  results  without 
re<iulrlng  inordinate  preparation  and  cooputer  time.  Although  the 
problem  has  not  been  resolved  in  this  testing  program,  a  possible 
rule  of  «•*»««■>«  for  ungaged  basins  is  to  model  all  the  streams  one 
order  less  than  the  main  stream  where  stream  order  is  determined 
by  the  pattern  of  confluences  of  tributary  streams  and  increases 
in  tile  downstream  direction.  This  will  probably  provide  a 
reasonable  representation  of  basin  slopes.  On  gaged  basins, 
lower  resolution  models  can  be  used  and  calibration  parameters 
can  be  adjusted  to  coiqiensate. 

S.2.1  Potter  Valiev  Watershed  Model 

As  with  the  Castro  Valley  watershed,  CH  losses  were  used  to 
model  rainfall  excess  in  the  initial  runs.  Since  Potter  Valley 
is  a  non-urban  basin,  Clark  Unit  Graph  method  and  Muskingum 
channel  routing  were  used  to  model  the  sub-basin  runoff  and 
channel  flow,  respectively.  The  HEC-1  input  data  for  the 
uncalibrated  Potter  Valley  model  are  on  file  at  the  HEC.  Figure 
21  shows  the  computed  and  observed  l^drographs  for  the  December 
20,  1964  flood  event  generated  using  the  uncalibrated  model. 

As  with  the  Castro  Valley  model,  the  timing  is  quite  a  bit 
slower  for  the  computed  hydrograph.  The  causes  of  this  slow 
response  are  different  for  the  Potter  Valley  model  however.  In 
the  Castro  Valley  model,  the  timing  problems  are  probably  a 
result  of  the  Inability  of  the  present  model  to  adequately  handle 
man-made  drainage  structures.  For  Potter  Valley,  the  problems 
appear  to  result  from  the  otethods  used  to  define  model 


coefficients . 
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FIGURE  21  POTTER  VALLEY  HECl -ADAPT  CLARK  MODEL 
UNCALIBRATED  DECEMBER  20,  1964  EVENT 


As  sencloned  previously,  the  Potter  Valley  DIM  did  not 
capture  all  tributaries  of  the  east  Fork  Russian  River  with  the 
result  that  aodel  slopes  are  less  than  actual  slopes  in  suse 
areas.  To  compound  this  problem,  the  overlaxid  flow  paths  derived 
by  the  model  can  be  quite  contorted.  (This  effect  is  explained 
in  a  following  paragraph.)  Another  source  of  error  is  the  way 
overland  and  channel  velocities  are  computed.  The  interface 
program  HECAD  uses  a  simplified  Manning's  equation  which  computes 
velocities  assuming  a  constant  depth  of  flow. 

The  overall  Impact  of  these  problems  is  seen  in  Figure  21. 
The  computed  hydrograph  peaks  seven  hours  later  than  the  observed 
hydrograph  and  the  peak  is  about  18  percent  less  than  the 
observed.  Runoff  volumes  are  similar.  Table  8  tabulates  these 
results.  Isolation  of  these  errors  to  determine  their  individual 
effect  on  the  simulation  results  is  discussed  below. 

In  order  to  gage  the  impact  of  topographic  model  errors  on 
the  hydrograph  timing,  HEC-1  input  data  were  derived  for  a 
Kinematic  Wave  model  using  the  same  loss  function  and  roughness 
as  for  the  Clark  BK>del.  Figure  22  shows  the  hydrographs  computed 
using  both  models.  The  spikes  of  the  Kinesuitic  Wave  hydrograph 
coincide  quite  well  with  the  spikes  in  the  observed  hydrograph 
leading  to  the  conclusion  that  the  timing  errors  are  mainly  the 
result  of  problems  with  the  derived  Clark  and  Muskingum 
coefficients  and  not  the  result  of  errors  in  the  topographic 
representation.  Possible  reasons  for  this  are  discussed  below. 

The  interface  program  HECAD  develops  time-area  c\irves  by 
detenslning  the  travel  time  of  each  triangular  element  in  the 


TABLE  8 


OOHPARISON  OF  COMPUTED  AND  OBSERVED  HYDROGRAPHS 
POTTER  VALLEY  HECl-ADAPT  CLARK  MODEL 
UHCAUBRATED  DECQIBER  20.  1964  EVENT 


SUM  OF 
FLOHS 

EQUIV. 

DEPTH 

(cfs-hr) 

(in) 

ConpuCed 

Hydrograph 

516747 

8.606 

Observed 

Hydrograph 

539922 

8.992 

DIFFERENCE 

-23175 

-0.386 

PERCENT 

DIFFERENCE 

-4.29 

-4.29 

MEAN 

TIME  TO 

PEAK 

TIME  TO 

FLOW 

CENTER 

OF  MASS 

FLOW 

PEAK 

(cfs) 

(hrs) 

(cfs) 

(hrs) 

6459 

49.63 

15320 

48.00 

6749 

40.51 

18700 

41-.  00 

-290  9.12 

-3380 

7.00 

-4.29  22.51 

-18.07 

17.07 

sub-basin  Co  tha  do«msCraaa  vertex  of  the  screen  link.  The  £io*i 

path  for  each  triangular  element  is  deterained  through  a 

centroid-to-centroid  routing  technique.  As  nentioned  previously 

this  routing  can  generate  flow  paths  that  are  quite  contorted  - 

with  s<MM  flow  paths  such  longer  chan  the  actual  flow  paths  (see 

Figure  23).  This  alone  can  cause  excessively  long  travel  tines. 

Additionally,  the  overland  and  channel  flow  velocities  of  each 

triangular  elenent  are  calctilated  using  a  sinplified  application 

of  the  Hanning's  equation  in  which  constant  depth  is  assuand. 

•  . 

This  can  cause  errors  in  Che  veloci^  conputation.  Since  the 
Clark  coefficients,  Tc  and  R,  and  the  Muskingun  coefficient,  K, 
are  also  derived  using  triangle  travel  tians,  the  errors  are 
conpounded. 

To  conpensate  for  these  errors,  calibration  of  the  Clark 
nodal  required  the  use  of  ssmII  roughness  factors  and  the 
inicial/unifom  loss  function.  Three  scorn  events  were  used. 
Results  of  this  calibration  are  shown  on  Figures  24,  25  and  26. 
HEC-1  input  data  for  the  calibrated  aradel  are  on  file  at  the  HEC. 
For  coaq>arison,  hydrographs  generated  by  HEC-1  using  optinlzed 
Clark  unit  graph  paraneters  are  also  shown  (17).  A  tabulated 
conparlson  of  the  observed  and  calibrated  laodel  hydrographs  for 
Che  Deceaber  20,  1964  stom  appears  on  Table  9. 

One  can  see  from  Figures  24  ,  25  and  26  that  the  HECl -ADAPT 
and  optinlzed  Clark  models  produce  hydrograph  peaks,  volumes  and 
tining  that  are  quite  sinllar  for  each  of  these  events.  After 
conpleting  the  calibration,  both  models  are  verified  using  the 
Jataiary  14,  1974  flood  event. 


Source:  R.N.  Eli,  "Runoff  and  Erosion  Predictions  Using  a  Surface 
Mine  Digital  Terrain  Model 


FIGURE  23  CENTROID-TO-CENTROIO  ROUTING 
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T  CLARK  FIGURE  24  POTTER  VALLEY  HECl -ADAPT  CLARK  MODEL 

CALIBRATED  (HEC-1  OPTIMIZED  CLARK  MODEL 
ALSO  SHOWN)  DECEMBER  20,  1964  EVENT 


00091 
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POTTER  HECl -ADAPT  CLARK  FIGURE  25  POTTER  VALLEY  HECl -ADAPT  CLARK  MODEL 

POHER  OBS  FLOW  .  CALIBRATED  (HEC-1  OPTIMIZED  CLARK  MODEL 

POTTER  OPT  CLARK  ALSO  SHOWN)  DECEMBER  18,  1955  EVENT 


TABLE  9 


COMPARIS(»l  OF  OOMPOTED  AMD  OBSERVED  HYDROGRAFHS 
POTTER  VALLEY  HECl -ADAPT  CLARK  MODEL 
CALIBRATED  DECEMBER  20.  1964  EV^ 


SUM  OF 

EQUIV. 

MEAN 

TIME  TO 

PEAK 

TIME  TO 

FIOUS 

DEPTH 

FLOW 

CENTER 

OF  MASS 

FLOW 

PEAK 

(cfs-hr) 

(In) 

(cfs) 

(hrs) 

(cfs) 

(hrs) 

ConpuCed 

Hydrograph 

522499 

8.702 

6531 

42.64 

21724 

43.00 

Observed 

^fdrogreph 

539922 

8.992 

6749 

40.51 

18700 

41.00 

DIFFERENCE 

-17423 

-0.290 

-218 

2.14 

3024 

2.00 

PERCENT 

DIFFERENCE 

-3.23 

-3.23 

-3.23 

5.27 

16.17 

4.88 

76 


Figure  27  shows  the  hydrographs  for  these  slonilations  %fhile 
Table  10  eoapares  Che  slsnilation  results.  Both  the  optialzed 
Clark  and  HECl-AOAPT  aodels  do  a  good  Job  of  reproducing  this 
event.  As  with  Che  calibration  events,  the  two  nodels  produce 
hydrographs  Chat  are  quite  similar.  From  Table  10  it  is  seen 
that  simulated  hydrograph  parameters  vary  by  about  13  percent. 

As  it  did  for  Castro  Valley,  Che  HECl-AOAPT  system  simulates 
rainfall -runoff  on  Potter  Valley  about  as  well  as  Che  HEC-1 
optimized  Clark  method. 

5.2.2  Sensitivity  Analysis 

The  results  of  Che  Potter  Valley  sensitivity  analysis  are 
quite  different  from  those  for  Castro  Valley.  For  the  Potter 
Valley  model,  peak  discharge  is  most  sensitive  to  overland 
roug^iness  and  channel  rouj^ess  parameters  while  lag  has  about 
the  same  sensitivity  to  all  calibration  parameters.  Equivalent 
depth  is  most  sensitive  to  the  estimates  of  loss  parameters. 

Table  11  CabxilaCes  these  results. 

The  differences  in  these  sensitivity  analyses  are  mainly  the 
result  of  variation  between  the  physical  characteristics  of  Che 
two  basins.  Runoff  does  not  concentrate  as  rapidly  in. Potter 
Valley  as  it  did  in  Castro  Valley  because  the  basin  is  ,much 
larger  and  non-urban.  Consequently,  the  estimation  of  the  loss 
function  parameter  becomes  less  significant  for  the  peak  while 
the  overland  and  channel  roughness  parameter  estimates  become 
more  significant. 

Channel  roughness  is  more  significant  in  the  Potter  Valley 
BKKiel  because  the  proportion  of  total  travel  time  accounted  for 


iseo* 
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TABLE  10 

C0MPARIS«»1  OP  POTTER  VALLEY  HECl-ADAFT  CLARK  HODEL 
AND  HEC-1  OPTIMIZED  CLARK  MODEL  HYDROGRAniS 
DECEMBER  20,  1964  EVENT 


SUM  OF 

BQOIV. 

MEAN 

TIME  TO 

PEAK 

TIME  TO 

FLOWS 

DEPTH 

FLOW 

CENTER 

OF  MASS 

FLOW 

PEAK 

(cfs-hr) 

(in) 

(cfs) 

(hrs) 

(cfs) 

(hrs) 

HBCl -ADAPT 

264064 

4.398 

3301 

45.07 

11914 

50.00 

Model 

Optimized* 

302303 

5.081 

3779 

44.53 

12574 

50.00 

Clark  Model 

.. 

DIFFERENCE 

-38239 

-0.683 

-478 

0.54 

-660 

0.0 

PERCENT 

•12.65 

-13.44+ 

-12.65 

1.21 

-5.25 

0.0 

DIFFERENCE 

• 

- 

OBSERVED 

307073 

5.161 

3838 

44.76 

11900 

50. ( 

*  Source:  HEC,  1984,  Spillway  Adequacy  Study  -  Coyote  Dam  and 
Lake  Mendocino. 

+  Discrepancy  between  sum  of  flows  and  equivalent  depth  Is  caused 
by  slight  differences  in  drainage  area  between  the  two  BK>dels. 


TABLE  11 


SEMSlTIVm  AMALYSIS 
POTTER  VALLEY  1R>DEL 


PARAMETER 

CASE 

PEAK 

LAG 

DISCH 

(%) 

(cfs) 

<hrs) 

Overland 

Ref. 

21724 

42.64 

Roughness 

•f20 

21030 

42.91 

+10 

21241 

42.81 

-10 

22217 

42.52 

-20 

22730 

42.36 

Channel 

Ref. 

21724 

42.64 

Roui^bneas 

+20 

20788 

43.06 

• 

+10 

21297 

42.84 

-10 

21938 

42.44 

-20 

22690 

42.23 

Initial/ 

Ref. 

21724 

42.64 

Uniform 

+20 

20948 

43.06 

Loss 

+10 

21337 

42.84 

-10 

22109 

42.43 

-20 

22501 

42.21 

BQUIV 

SENSITIVITY  COEFFICIENTS 

DEPTH 

a9Z9* 

al/l* 

AD/D* 

AB/B 

AB/B 

AB/B 

(in) 

8.702 

8.631 

0.160 

0.032 

0.041 

8.666 

0.222 

0.040 

0.041 

8.742 

0.227 

0.028 

0.046 

8.786 

0.232 

0.033 

-0-.048 

8.702 

8.685 

0.215 

0.049 

0.010 

8.694 

0.197 

0.047 

0.009 

8.713 

0.099 

0.047 

0.013 

8.721 

0.222 

0.048 

0.011 

8.702 

8.080 

0.179 

0.049 

0.357 

8.388 

0.178 

0.047 

0.361 

9.038 

0.177 

0.049 

0.386 

9.931 

0.179 

0.050 

0.396 

*  jVQ  *”  change  in  peak  discharge 

Q  -  reference  value  of  peak  discharge 
AL  -  change  in  lag 
L  —  reference  value  of  lag 
AD  *  change  In  e<tulvalenC  depth- 
D  -  reference  valiM  of  equivalent  depth 
AB  -  change  in  paraneter  (e.g.  overland  roughness) 
B  —  reference  value  of  parameter 
AQ/Q  ••  sensitivity  coefficient  for  peak  discharge 
AB/B 


Hote;  reference  values  for  state  variables  and  parameters 
are  from  the  calibrated  model 


by  ehaniMl  flow  is  sure  on  larger  basins  than  on  saaller  basins. 
Potter  Valley  also  has  fairly  steep  overland  topography  resulting 
in  shor.'ei*  overland  lag  tines  relative  to  the  less  steep  channel 
segnents . 

The  proportioning  effect  is  also  seen  in  the  sensitivity  of 
the  lag.  For  Castro  Valley,  the  lag  was  nost  sensitive  to 
overland  roughness,  while  in  the  Potter  Valley  model,  lag  has 
about  the  sane  sensitivity  for  both  overland  and  channel 
roughness . 

As  with  the  Castro  Valley,  equivalent  depth  is  most 
sensitive  to  the  loss  function  parameter  estimate  which 
determines  the  volume  of  runoff. 

S.3  Program  Modifications  to  Improve  Results 

Modifications  to  Improve  the  modeling  of  urban  basins  are 
suggested  in  the  recommendations  section  of  this  report.  A 
simple  modification  to  Improve  the  modeling  of  natural  basins  is 
tested  and  the  results  discxissed  below. 

The  effects  of  basin  and  channel  storage  are  important  in 
natural  basins.  In  HECl-ADAPT,  the  options  available  to  model 
storage  effects  are  the  SCS  Dimensionless  Unit  Graph,  the  Clark 
and  the  Snyder  Unit  Graphs  and  the  Muskingum  and  normal>depth 
methods  for  channel  routing.  The  coefficients  for  each  of  these 
options  (except  the  normal-depth  option)  are  derived  through  the 
computation  of  triangle  travel  times.  The  problems  associated 
with  the  triangle  travel  time  computation  were  discussed 
previously.  In  HECAD,  the  computation  of  overland  velocities  is 
made  in  one  statement  of  the  overland  routing  subroutine. 


Modificaeion  of  this  sxibroutina  was  tharafora  8laq>la  and  is 
discussad  balow. 

The  SCS  (30)  has  davalopad  a  chart  ralaCing  land  usa,  slopa 
and  overland  valoclty  (saa  Flgura  28).  In  a  Cast  aodlflcation, 
this  infomation  has  haan  convartad  to  aquation  fora  and 
incorporated  into  HECAD.  The  HEC-1  input  data  derived  using  tiiis 
version  of  HECAO  are  on  file  at  the  HEC.  Figure  29  shows  the 
hydrographs  generated  by  the  uncalibrated  Potter  Valley  model 
with  and  without  the  new  overland  valpci^  routine. 

^drograph  timing  is  anich  closer  to  tlie  observed  for  die 
hydrograph  generated  using  the  modified  routine.  The  quicker 
response  has  also  concentrated  more  of  the  runoff  volume  within 
the  80  minute  time  base  maiding  both  the  peak  discharge  and 
equivalent  depths  greater  than  those  of  the  previously  computed 
hydrograph.  Although  the  results  are  not  conclusive,  this  simple 
modification  appears  to  yield  a  sig^ficantly  better  response.  A 
similar  modification  may  be  possible  for  the  computation  of 
stream  velocities  since  the  Muskingum  coefficient  K  and  the 
number  of  time  steps  are  also  derived  using  a  simplified 
Manning's  equation  approach. 


6.  Reconaendations 


TWO  sets  of  recooBsendacions  are  given  in  this  section.  The 
first  suggests  changes  in  or  additions  to  ADAPT  software  i^ich  - 
will  nake  GIS  development  easier.  The  second  suggests  changes  or 
additions  to  HECAO  software  to  improve  sradeling  results. 

6.1  ADAPT  Changes .  The  HECl-ADAPT  system  consists  of  about 
forty  prograsus  and  a  library  of  utility  subroutines.  Of  these 
for^,  less  than  thir^  were  actually  used  in  the  testing  program 
and  many  of  these  programs  were  quite  small.  To  make  software 
management  easier,  this  set  of  programs  could  probably  be 
combined  into  a  single  program.  R.D.  Carl,  a  Hydraulic  Engineer 
in  the  Planning  and  Analysis  Branch  of  the  HEC,  performed  a 
preliminary  assessment  of  ADAPT  (3).  In  this  assessment,  he 
suggests  that  a  main  program  be  developed  to  connect  the  programs 
and  manage  the  user's  input  and  output.  The  management  routines 
could  contain  additional  error  checking  and  data  validation 
procedures  and  provide  guidance  to  the  user  to  help  him  follow 
the  flow  diagrams  contained  in  the  ADAPT  documentation. 

A  major  problem  with  the  existing  sofware  is  that  no 
consistent  format  is  used  for  entering  data  to  the  various 
programs.  In  his  assessment,  R.D.  Carl  also  suggests 
incorporating  a  format  similar  to  that  used  by  other  HEC  programs 
which  employ  a  record  Identifier  (3).  This  single  improvement 
would  speed  GIS  development  considerably. 

An  Important  part  of  TIN  error  correction  is  accomplished 
using  plots  of  the  base  map,  stream  networks,  triangle  slope 


dlr«ccions  atid  bastn  contours.  A  hand*autde  plot  of  the  overland 
flow  network  has  also  been  found  to  be  quite  useful  for  error 
correction.  This  plot  is  presently  constructed  ixsing  output  fron 
the  ADAPT  prograa  LISNET.  Addition  of  this  plotting  capability 
would  speed  the  error  correction  process. 

The  digitizing  software  and  hardware  are  very  ioiportant  to 
the  development  of  a  CIS.  During  this  testing  program,  the 
normal  digitizer  setup  had  to  be  altered  to  extable  the  use  of  the 
computer  program  DIGITZ  (which  writes  digitizer  output  ixn:o*a 
file  usltig  the  proper  format)  tdxile  digitiziitg.  This  hardware 
problem  needs  to  be  corrected.  The  digitizing  software  should 
be  Improved  to  sudee  it  easier  to  input  soil  and  laxtd  use 
polygons.  At  present,  the  triangle  input  software  is  xised  to 
generate  input  data  describing  these  polygons.  The  input  formats 
required  by  the  varloxis  programs  are  different.  Consequently, 
output  from  the  digitizing  must  be  hand  edited. 

Development  of  a  GIS  requires  assigning  soil  and  land  use 
to  each  triaxigle.  In  oMSt  cases  it  is  much  easier  to  hai\dle  the 
distribution  of  soil  and  land  use  types  in  the  GIS  by  assigning  a 
percentage  mix  to  each  triangle.  If  this  is  not  done,  the  DTK 
must  be  modified  to  Incorporate  the  soil  and  land  use  boundaries. 
This  can  Involve  considerable  work  redefining  triangles  and 
making  sure  that  the  topographical  representation  remains  intact. 
The  present  model  allows  assignment  of  a  mix  of  land  uses  to 
individual  triangles  through  the  program  LDIN.  This  same 
capability  needs  to  be  developed  for  inserting  mixes  of  soil 


types  by  triangle. 


6.2  HECAD  Changes .  Given  an  established  CIS,  there  are  nany 
ways  to  derive  input  for  hydrologic  aodeling.  In  HECl -ADAPT  for 
exaiqtle,  CN's  are  derived  from  associations  of  land  use  and  soil 
type  in  each  triangle.  Clark  coefficients  are  derived  using 
triangle  slopes,  roughnesses  and  slope  directions.  Alternate 
methods  for  computing  these  coefficients,  like  those  described  by 
Ell  (6,7)  and  Li  (25),  are  available  and  could  be  incorporated 
into  tlie  existing  program  to  take  better  advantage  of  the 
capabilities  of  the  CIS. 

Eli  has  developed  a  methodology  that  computes  routing 
coefficients  using  the  overland  flow  paths  defined  by  triangle 
slopes  (6).  Incorporation  of  this  method  or  a  similar  method 
would  probably  improve  the  derivation  of  Clark,  Snyder,  SCS  and 
Muskingum  coefficients. 

An  alternate  method  of  computing  overland  velocities  was 
examined  earlier  in  this  report.  The  preliminary  testing 
suggested  that  this  method  may  be  better  than  the  one  presently 
employed. 

Only  two  methods  to  compute  infiltration  are  available  in 
the  present  model.  It  would  be  fairly  easy  to  add  a  Holtan 
method  option.  Additional  data  required  by  this  method  could  be 
stored  in  the  Soil  Matrix  File.  An  additional  HECAD  routine 
would  have  to  be  developed  to  derive  the  Holtan  parameters  and 
output  the  HEC-1  input  data  file. 

An  Improvement  in  the  modeling  of  urban  basins  would 
probably  be  accomplished  if  all  the  flow  elements  allowed  in  the 


HBC-l  Kineaatlc  Wave  option  were  utilized,  ie:  two  overland  flow 
el^aenta  aM  three  channel  elements.  The  additional  information 
required  could  be  stored  in  the  triangle  file.  Overland  flow 
elements  could  model  runoff  from  different  land  uses  in  each  sub 
basin.  Channel  elements  could  model  the  local  drainage 
systems,  interceptors  and  main  channels.  The  present  model 
overland  flow  plane  and  one  main  channel. 


allows  only  one 


7.  Suaaary  and  Conclusions 


7.1  Model  Resolution.  Although  only  United  testing  was 
done  on  the  effect  of  model  resolution,  the  results  appear  to  - 
conflm  those  of  earlier  researchers  (8).  Hl^er  resolution 
models  capture  more  of  the  existing  terrain  features  Including 
the  smaller  tributary  streams  and  valleys  and  provide  a  more 
accurate  topographic  -representation.  The  result  of  this 
Increased  resolution  Is  a  quicker  runoff  response  and  a  greater 
peak  discharge  on  the  basins  tested. 

As  stated  In  the  Potter  Valley  section  of  this  report,  a 
major  difficulty  for  the  user  Is  deciding  what  model  resolution 
Is  appropriate  for  the  purposes  of  a  given  study.  The  simple 
rule  of  thxinb  proposed  In  this  report  will  probably  be  more  than 
adequate  for  ralnfall^runoff  modeling  on  larger  basins. 

Restated,  the  rule  Is  to  model  all  streams  one  order  less  than 
the  main  stream  where  stream  order  Increases  In  the  downstream 
direction.  Additional  experience  with  this  technique  will  be 
necessary  before  a  more  definitive  solution  is  fotind. 

7.2  Modeling  Ungaged  Basins.  The  tuncallbrated  Castro 
Valley  and  Potter  Valley  models  produced  hydrographs  that  were 
quite  similar  to  the  historical  events.  In  both  situations  the 
predicted  volumes  were  within  about  six  percent  of  observed 
volumes.  Peak  discharges  were  under 'predicted  about  six  percent 
In  the  Castro  Valley  model  and  by  about  18  percent  In  the  Potter 
Valley  model.  In  both  models,  the  lag  times  were  over  20  percent 
greater  for  the  predicted  hydrographs  than  for  the  observed 


hydrographs.  This  disparity  in  lag  tine  is  the  najor  difference 
between  computed  and  observed  hydrographs  in  both  simulations. 

In  the  Castro  Valley  model,  the  timing  problems  are 
probably  due  to  the  Inability  of  the  model  to  account  for  the 
effect  of  man-made  drainage  structures  on  the  runoff  response. 
This  is  compensated  for  in  the  calibration  by  usin^  low  overland 
and  channel  •  roughness  factors  and  by  using  the  initial/uniform 
loss  function  to  get  more  volume  on  the  rising  side  of  the 
hydrograph.  In  the  Potter  Valley  model,  the  long  lag  is  .probably 
due  to  problems  with  the  way  Clark  and  Muskingum  coefficients  and 
the  time-area  curves  are  derived.  As  with  the  Castro  Valley 
model,  modified  roughness  factors  and  the  Initial/uniform  loss 
function  are  \tsed  to  obtain  a  calibration. 

Even  with  these  problems  the  simulation  results  are 
reasonable  for  ungaged  basins.  This  is  very  encouraging  because 
all  input  data  for  the  runoff  model  is  derived  directly  from  non- 
calibrated  CIS's  of  the  drainage  basins.  There  appears  to  be 
great  potential  for  CIS's  to  provide  a  basis  for  modeling  ungaged 
basins. 

7.3  Modeling  Urban  and  Non-Urban  Basins.  The  flexibility 
of  the  system  is  demonstrated  by  the  relative  success  in 
reproducing  hydrographs  for  both  an  urban  basin  and  a  non-urban 
basin.  Although  some  problems  remain,  in  general  the  methodology 
used  to  develop  coefficients  for  the  rainfall-runoff  model 
appears  to  -be  sound. 

In  conclusion,  the  HECl -ADAPT  system  accounts  for  the 


hydrologic  dlversl^  of  «  drainage  basin  and  accomodates  the 
derivation  of  runoff  model  coefficients.  Since  this  derivation 
is  based  on  a  physical  representation  of  the  basin,  HECl -ADAPT 
provides  a  reasonable  SMthod  for  simulating  the  response  of 
ungaged  watersheds.  The  methodologies  used  to  develop  both  the 
CIS  and  the  input  data  for  the  rainfall- runoff  model  are  fairly 
sound.  PrelisdLnary  analysis  and  experience  with  the  model 
indicate  that  slisple  modifications  to  both  ADAPT  and  to  the 
Interface  HECAD  could  be  implemented  and  would  improve  the  model 


results. 
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APPENDIX  A 


Modifications  and/or  Corrections  to  ADAPT  System  Routines 

The  first  two  tasks  of  the  testing  program  were  to  determine 
whether  the  necessary  hardware  and  software  to  develop  a  data 
base  were  available  and  whether  everything  was  working  properly. 

It  was  felt  that  this  could  be  best  accomplished  by  building  a 
CIS  from  scratch.  A  small  urban  basin,  Castro  Valley,  was  chosen 
for  this  purpose.  Castro  Valley  has  been  used  in  a  number  of 
other  studies  by  the  NEC. 

The  first  Castro  Valley  model  was  constructed  by  following 
Figure  1  in  Section  2,  Volume  1  of  the  HEC1>ADAPT  doctimentation. 
Program  errors  and  problems  were  corrected  as  they  were 
encoxmtered.  Some  bugs  may  remain  in  the  HECl -ADAPT  prograais  that 
were  not  used  in  this  test  application. 

The  HECl -ADAPT  programs  were  compiled  and  linked  by  R.D. 

Carl  at  the  HEC  in  January,  1985.  One  additional  program, 

DIGITZ,  was  compiled  and  linked  in  July,  1985.  The  DIGITZ 
program  provides  a  systematic  method  for  inserting  triangles 
using  the  HEC  digitizer.  A  short  description  of  this  program  and 
the  other  programs  used  in  the  testing  program  is  given  in 

f 

Table  1. 

As  mentioned  previously,  program  errors  and  problems  were 
corrected  as  required  during  the  process  of  building  a  GIS  for 
the  Castro  Valley.  Table  2  lists  these  program  changes.  A  hard 
copy  of  tiie  source  code  of  each  of  these  programs  is  on  file  at 
the  HEC  with  the  required  changes  to  the  original  code  indicated. 


After  coflq>leting  the  above  tasks  it  was  necessary  to  check 
the  cooputation  routines  within  the  interface  prograai  HECAD. 

Much  of  this  work  involved  inserting  %rrite  statements  into 
various  subroutines  of  the  interface  to  print  out  intermediate 
and  final  results.  These  results  were  then  verified  by  hand 
computations.  An  error  in  the  overland  travel  time  computation 
was  discovered  and  corrected  in  this  step. 

Other  subroutines  were  checked  by  comparing  HECAD  output 

with  Information  contained  in  the  data  base  and  listed  for  this 

« 

purpose.  The  HECAD-derived  Initial/uniform  loss  coefficients  and 
curve  numbers  were  checked  this  way.  Table  3  gives  a  list  of  the 
HECAD  routines  that  were  checked  and  modified.  The  modifications 
are  described  in  Table  2. 

Miscellaneoxis  Quirks  and  Errors 

The  following  section  describes  other  problems  or  possible 
errors  in  logic  that  were  encountered  while  developing  the  CIS 
for  this  testing  program. 

1.  There  is  a  problem  with  simulating  small  events  in  which 
the  initial  abstraction  is  greater  than  total  rainfall  on  some 
sub-basins.  HEC-1  will  not  run  when  Initial  abstraction  is 
greater  than  total  rainfall.  Thus,  HEC-1  will  abort  upbn 
encountering  a  sub-basin  with  less  rainfall  than  initial 
abstraction.  With  the  existing  setup,  flow  cannot  be  routed 
throu^  a  non-rainfall-excess  producing  basin. 

2.  The  horizontal  distance  rather  than  slope  distance  is 
used  for  computing  travel  times  for  the  Clark,  SCS,  and  Snyder 


■echods.  This  will  produce  errors  if  steep  ground  slopes  exist  in 
the  basin  being  studied. 

3.  HECAD  does  not  make  use  of  the  land  use  information  that 
is  established  in  the  triangle  file  using  the  programs  POLDIG, 
POLYCR,  and  UHPIN.  The  data  Inserted  using  these  programs  can  be 
accessed  by  FENPLT  to  make  Unique  Attribute  Polygon  (UAP)  and 
Single  Attribute  Polygon  (SAP)  plots,  but  it  is  still  necessary 
to  nm  program  LXJIN  to  insert  the  land  use  data  required  to  run 
HECAD. 

POLDIG,  POLYCR  and  UHPIN  can  also  be  bypassed  for  insertion 
of  soils  data.  Programs  ADSOIL,  which  establishes  the  percent  of 
each  hydrologic  group  associated  with  each  soil  type,  and  TRIINl, 
which  can  be  used  to  Insert  a  soil  number  into  a  given  column  of 
the  triangle  file,  can  be  run  to  write  this  information  into  the 
proper  files. 

Programs  POLYCR,  POLDIG  and  UNPIN  need  to  be  x:\in  only  when 
UAP  or  SAP  plots  are  desired. 

4.  In  the  triangle  file,  HECAD  requires  land  use  data  in 
Coliimns  30  through  33  and  a  soil  number  in  Column  49  in  order  to 
execute.  Tne  programs  described  in  Item  3  above  are  tised  to 
establish  this  data.  Additionally,  HECAD  requires  the  Hanning's 
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roughness  for  each  stream  link.  Program  NETIN  is  used  to  input 
the  roughness  into  Colximn  33  of  the  network  file.  None  of  these 
requirements  are  mentioned  in  the  HECl- ADAPT  documentation. 

5.  There  appear  to  be  some  problems  with  the  cuirve  number 
adjustment  algorithm.  For  antecedent  precipitation  values  lower 
than  about  0.9  inch,  the  adjusted  curve  number  is  actually  lower 


than  the  unadjusted  value  for  the  dormant  season.  This  doesn't 
seem  to  make  sense.  Otherwise,  the  algorithm  appears  to  give 
reasonable  values. 

6.  The  time -step  size  chosen  for  the  computation  interval 
should  be  a  function  of  the  stream  link  with  minimum  travel  time. 
If  the  computation  interval  is  larger  than  the  time  it  takes  for 
water  to  traverse  the  length  of  the  stream  link,  then  some  error 
will  be  introduced  into  the  routing  computations. 

Time  Required  to  Develop  a  CIS 

\ 

The  tiob  required  to  develop  a  GIS  is  of  course  a  function 
of  the  size  and  complexity  of  the  drainage  basin.  To  provide 
some  idea  of  this  time  requirement,  a  time  log  was  kept  during 
the  development  of  the  Castro  Valley  82- triangle  model.  Table  4 
shows  this  time  log. 

This  particular  data  base  took  approximately  one  week  to 
develop.  However,  a  major  requirement  not  included  in  this  log 
is  the  time  to  get  "up  to  speed"  with  the  technology.  The  first 
Castro  Valley  model  provided  a  simple  data  base  for  learning  the 
structure  and  procedures  of  the  method.  Thus,  the  one-week 
period  to  build  the  82 -triangle  model  assumes  prior  experience 
with  HECl -ADAPT. 

Documentation 

Table  5  tabulates  the  HECl-ADAPT  programs  used  in  the 
testing  program,  the  source-  and  executable  file  names  and  the 
input  and  output  files  required  to  run  them  for  the  82 -triangle 
Castro  Valley  model. 
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The  computer  code  sheets  that  were  used  to  run  the  programs 


are  also  on  file  at  the  HEC. 


TABLE  1 


PROGRAM 

DIGITZ 

DIGIT 

FILEST 

VERTEX 

CRETRI 

PENPLT 

FIXVER 

FIXTRI 

BOUNV 

ADJCHK 

OtNET 

EDITNT 

OVERLN 

LISHET 

VUST 

TLIST 

POLDIG 

UNPIH 

POLYCR 

IMTSMF 

ADSOIL 


PROGRAM  DESCRIPTIONS 


DESCRIPTION 


Provides  a  sysCematlc  way  of  inputting  triangle  data 
tising  the  HEC  DATATB  II  digitizer. 

Prepares  a  triangle  data  file  from  the  digitizer  data 
file  for  input  to  FILEST. 

Converts  output  of  DIGIT  to  input  for  VERTEX  and 
CRETRI. 

Establishes  vertex  file  using  input  generated  by 
FILEST  or  nanually. 

Establishes  triangle  file  and  adds  topology  data  to 
vertex  file. 

Produces  display  plots  of.  GIS  at  any  scale.  Also 
used  for  error  checking. 

Corrects  or  aK>difies  topology  in  data  base  by 
deletion,  redefinition  or  addition  of  vertices. 
Corrects  or  nodifies  topology  in  data  base  by 
deletion,  redefinition  or  addition  of  triangles. 
Establishes  boundary  file  which  is  input  for  EDITNT 
and  OVERIH  programs.  Also  used  to  identify  topologic 
hole  problems. 

Identifies  adjacency  and  topologic  errors  in  the 
data  base. 

Establishes  stream  network  and  stream  drainage  file. 

Corrects  topologic  errors  in  stream  file.  Also 
Identifies  topographic  errors  in  data  base. 
Establishes  overland  drainage  network  and  stores 
results  in  stream  network  file. 

Produces  a  listing  of  the  values  of  selected  data 
types  stored  within  the  stream  network  file. 

Produces  a  listing  of  the  vertex  file. 

Produces  a  listing  of  the  triangle  file. 

Inserts  UAP  numbers  or  SAP  values  in  a  specified 
column  of  the  triangle  file  using  digitizer 
coordinates. 

Assigns  attribute  values  to  triangles  based  on  the 
UAP  nunbers  they  have  been  assigned. 

Allows  manual  definition  of  a  vertex  chain  for 
Insertion  of  SAP  or  UAP  values.  Can  also  be  \ised  for 
coorecting  POLDIG  errors. 

Initializes  the  soil  matrix  file. 

Used  to  modify  soil  charateristic  values  in  the  soil 
matrix  file. 
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TABLE  I  (Cent.) 

PROGRAM  DESCRIPTION 

LOIN  Used  to  insert  the  percentage  mix  of  land  uses  in  each 

triangle  into  the  triangle  file. 

TRIIMI  Used  to  change  or  insert  a  nunber  or  value  in  a 

specified  colunn  of  the  triangle  file  (e.g.  used  to 
insert  soil  number  in  col.  49  in  this  testing 
program) . 

NETIN  Used  to  Insert  stream  link  attribute  values  into  the 

stream  network  file  (e.g.  used  to  insert  Manning's 
roughness  for  each  stream  link  into  col.  33  in  this 
testing  program) . 

SEGINS  Used  to  insert,  sub-watershed  identifiers  for  each 

stream  link  into  the  stream  network  file. 

HECAD  Interface  program  which  connects  ADAPT  and  HEC-1. 


HEC-1 


The  HEC  rainfall -runoff  program. 
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PROGRAM 

DIGIT 

POTL 

POLDIG 

PENPLT 

OVERLN 

HETIM 

SEGIMS 

HECAD 


TABLE  2 

PROGRAM  MODIFICATIONS 


MODIFICATION 

Modified  the  call  sCateisent  "CALL  STATNE" .  STATNE  is 
a  subroutine  in  the  library  of  stibroutines  which 
computes  easting  and  northing  based  on  a  given 
latitude,  longitude  and  zone.  The  call  statement 
specified  the  wrong  variable  name  for  zone. 

Modified  a  read  statement  in  FUTL.  PUTL  is  a. 
library  subroutine  that  is  called  during  the 
execution  of  CRENET.  CRENET  was  aborting  on  an  "end- 
of-file"  error.  Inserted  an  "End-  "  into  the  read 
statement. 

Modified  two  read  statements  to  allow  use  of  the 
digitizer  program  DIGITZ  for  insertion  of  polygon 
data. 

Modified  a  write  statement  in  subprogram  CONTOR.  The 
original  subprogram  specified  the  wrong  logical  file 
number. 

Added  a  common  block  to  each  of  these  programs.  The 
library  subroutine  GETL  is  called  during  the 
execution  of  each  program  and. requires  the  variables 
in  the  common  block. 

Changed  a  variable  name  in  subroutine  GLAND  to 
correct  an  error  in  ^e  travel  time  computation. 
Modified  GLAND  and  CHANK  subroutines  to  get  proper 
insertion  of  the  KG  card  in  the  HEC>1  data  file  for 
the  kinematic  wave  option. 

Modified  GLAND  to  compute  cumulative  time -area  curve 
coordinates.  Grlginal  subroutine  computed  incremental 
ordinates.  This  resulted  in  HEC-1  computing  negative 
unit  hydrograph  ordinates. 

Modified  GLAND  to  do  Interpolation  on  cumulative 
time -area  curve  ordinates  to  obtain  a  smoother 
function. 

Modified  GLAND  to  correct  travel  time  vs.  stream 
length  histogram.  Grlginal  computation  left  out  the 
first  two  links  in  the  non-klnematic  wave  runs. 
Modified  rain  gage  wei^t  computation  routines  in 
subroutine  RAINW.  Grlginal  routine  did  not  work. 
Modified  subroutine  GLAND  to  get  baseflow  cards  for 
Clark,  Snyder  and  SCS  methods.  Grlginal  routine 
inserted  only  baseflow  cards  for  the  kinematic  wave 
option. 


TABLE  3 


HECAD  PROGRAH  CHECKS 


HECAD  COMPUTATIONS 

CHECKED 

MODIFIED 

Rain  gage  weighting 

temporal 

* 

* 

areal 

* 

*  (new  routine) 

Loss  parameters 

Inltlal/unlfom 

* 

curve  nuiriiers 

* 

Curve  nuober  adjustment 

* 

Percent  Imperviousness 

Overland  flow  parasMters 

Clark 

* 

*  (time-area  curve) 

Snyder 

*  (time-area  curve) 

SCS 

* 

*  (time-area  curve) 

Kinematic  wave 

* 

*  (Insertion  of  KO 

card) 

Channel  routing  parameters 

Muskingum 

* 

Kinematic  wave 

* 

Baseflow  parameters 

* 

*  (Insertion  of 

baseflow  cards 
for  Clark,  Snyder 
and  SCS) 


Note:  unlfon  flow,  reservoir  touting  and  channel  loss 
coiq>utatlons  were  not  checked. 
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TABLE  4 

CHRONOLOGY  OF  THE  82 -TRIANGLE 
CASTRO  VALLEY  MODEL  CONSTRUCTION 


DATE 

TIME 

COMMENT 

September  20,  1985 

9:00-16:30 

Create  and  digitize  TIN. 

September  23 ,  1985 

9:30 

DIGIT 

m  mm 

10:30 

FILEST 

m  mm 

11:00-11:30 

VERTEX 

m  mm 

14:00 

CRETRI 

m  mm 

15:00 

PENPLT 

M  mm 

16:00-16:30 

FIXVER 

Septeidier  24,  1985 

8:00 

PENPLT 

m  mm 

10:00 

BOUNV 

m  mm 

10:15 

ADJCHK 

m  mm 

10:30 

CRNET 

m  mm 

11:00 

OVERIM 

m  mm 

11:30-12:30 

USNET 

September  25 ,  1985 

12:00-16:30 

TIN  modification.  Use  DIGIT 
to  determine  coordinates  of 
new  triangles. 

September  26,  1986 

8:00 

FIXVER 

m  mm 

9:00 

FIXTRI 

m  mm 

10:00-12:00 

BOUNV,  ADJCHK,  CRNET.  OVERIN 
LISNET 

m  mm 

14:00-15:30 

PENPLT 

m  mm 

15:30-17:00 

POLDIG  input  preparation 

Septeiid>er  27,  1985 

9:00-11:00 

POLDIG 

m  mm 

15:30-16:30 

POLYCR 

September  30,  1985 

9:00-10:00 

PENPLT 

m  mm 

10:30-11:30 

UNPIN,  PENPLT 

m  mm 

15:00-16:30 

INTSMF,  ADSOIL,  LUIN 

m  mm 

16:30-17:00 

TRIINl 

October  1,  1985 

8:00-10:00 

NETIN 

m  mm 

13:00-13:30 

SEGINS 

m  mm 

13:30-14:00 

HECAD 

m  mm 

14:00-15:00 

HEC-1 

Total 


about  36  hours 


TABLE  S 


adapt  programs  used  to  develop  82-triangle 

CASTRO  VALLEY  MODEL 


PROGRAM  SOURCE  EXECUTABLE 

DIGIT  DIGIT. S  DIGIT. X 


FILEST  FILEST.S  FILEST.X 


VERTEX  VERTEX.  S  VERTEX.X 


CRETRI  CRETRI.S  CREIRI.X 

PENPLT  PENPLT.S  PENPLT.X 


FIXVER  FIXVER.S  FIXVER.X 


DATA 

JOBSTREAM 

1  :  header  cards 

DIGIT. J 

TAPE20  (digitizer 

deck) 

o  :  DIGIT. 2 

OUTDIG  (printer 

output) 

i  :  DIGIT. 2 

FLEST.J 

FILEST. 2 

•  . 

o  :  FLT.2 

FLV.2 

FLP.2 

FILEOUT  (printer 

output) 

i  :  VRTEX.2 

VRTEX.J 

FLV.2 

FLP.2A  (FLP.2  plus 

header  card) 

o  ;  CS'-RO.Vl 

VERTOUT  (printer 

output) 

i  :  FLT.2A 

CRTRI.J 

o  :  CSTRO.Tl 

CSTRO.Vl 

CRTOUT  (printer 

output) 

i  :  PEN.l 

PEN.J 

CSTRO.Vl 

CSTRO.Tl 

o  :  PLT.P 

PENOUT  (printer 

output) 

i  :  FXVR.2 

EXVR.J 

o  :  CSTRO.Vl 

CSTRO.Tl 

FXVROUT  (printer) 

Note:  i  -  input  file 
o  -  output  file 


TABLE  5  <cont.) 


PROGRAM  SOURCE  EXECUTABLE  DATA 


BOUNV  BOUNV.S  BOUHV.X 


i  :  CSTRO.Tl 
o  :  CSTRO.Bl 

BOUNOUT  (printer} 


AOJCHK  ADJCHK.S  ADJCHK.X 


i  :  CSTRO.Tl 
CSTRO.Vl 

o  :  ACHKOUT  (printer) 


CRNET  CRNEI.S  CRNET.X  i  :  CSTRO.Tl 

CSTRO.Vl 
o  :  CSTRO.Nl 

CRNTOUT  (printer) 


USNET  LISNET.S  LISNET.X 


i  :  LSNT.U 
LSNT.2A 
CSTRO.Ml 
CSTRO.Tl 

o  :  LSNTOUT  (printer) 


OVERLN  OVERIH.S  OVERLN.X 


i  :  CSTRO.Bl 
CSTRO.Vl 
CSTRO.Tl 
o  :  CSTRO.Nl 

OVROUT  (printer) 


DIGIT  DIGIT.S  DIGIT.X 
Note:  used  second  tine  to 
deteraine  coordinates  of  new 
triangles  added  to  data  base 


1  :  header  cards 
TAFE20 
o  :  DIGIT. 3 

OUTDIG  (printer) 


FIXVER  FIXVER.S  FIXVER.X 


i  :  FXVR.3 
o  :  CSTRO.Vl 
CSTRO.Tl 

FXVROUT  (printer) 


FIXTRI 


FIXTRI.S  FDCTRI.X 


i  :  FXTR.2 
o  :  CSTRO.Vl 
CSTRO.Tl 

FXTROUT  (printer) 


BOUNV  * 

AOJCHK  * 

CRNET  same  as  before 
OVERLN  * 

USNET  * 


JOBSTREAM 
BOUNV. J 

AJCHK.J 

CRNET. J 

LSNT.J 

OVR.J 

DIGIT. J 

FXVR.J 

FXTR.J 


.'ABLE  S  (cont.) 


PROGRAM  SOURCE  EXECUTABLE 

POLDIG  POLDIG.S  POLOIG.X 


POLYCR  POLYCR.S  POLYCR.X 

UNPIN  UNPIN. S  UNPIN. X 


INTSMF  INTSMF.S  INTSMF.X 

ADSOIL  ADSOIL.S  ADSOIL.X 

LUIN  LUIN.S  LUIN.X 

TRIINl  TRIINl.S  TRIINl.X 

NETIN  NETIN.S  NETIN.X 

BEGINS  BEGINS. S  SEGINS.X 


data  JOBSTREAM 

i  :  header  card  PLDG.J 

TAPE20 
CSTRO.Vl 
CSTRO.Tl 
o  :  PLDG.2A 

PLDGOUT  (printer) 

i  :  PLCR.2  PLCR.J 

CSTRO.Tl 
CSTRO.Vl 

o  :  PLCROUT  (printer) 

i  :  UNPN.IA  UNPN.J 

UNPN.3A 
UNPN.4A 
CSTRO.Tl 
o  :  UNPN.2A 

UNPNOUT  (printer) 

i  :  INT.2  INT.J 

o  :  SMF.2 

INTOUT  (printer) 

i  :  AD,2  AD.J 

SMF.2 

o  :  ADOUT  (printer) 

i  :  LUIN.2  LUIN.J 

CSTRO.Tl 

o  :  LUINOUT  (printer) 

i  :  TRN.2  TRN.J 

CSTRO.Tl 

o  :  TRNOUT  (printer) 

i  :  NTN.2  NTN.J 

CSTRO.Nl 

o  :  NTNOUT  (printer) 

i  :  SGN.2  SGN.J 

CSTRO.Nl 
CSTRO.Tl 

o  :  SGNOUT  (printer) 


TABLE  S  (cone.) 


PROGRAM  SOURCE  EXECUTABLE  DATA 

HECAD  HECAD.S  HECAD.X  1  :  CAL.l 

RES.l 
DIV.l 
RAIN.l 
HCD.l 
SMF.2 
CSTRO.Nl 
CSTRO.Tl 
o  :  HCD.2 

HCDOUT  (printer) 

i  : 


JOBSTOEAM 

HCD.J 


HEC-1 


HECl 


o 


HCD.2 

header  cards 
HECIOUT  (printer) 


HECl.J 


